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ПАМЯТИ УЧЕНОГО 





21 июня 2020 года на 76-м году жизни скоропостижно скончался 
замечательный ученый, инженер и организатор, доктор технических наук, 
профессор Рудольф Анатольевич Нейдорф. 17 октября ему исполнилось бы 76 лет. 
Рудольф Анатольевич Нейдорф — известный в России и за рубежом специалист в 
области автоматических и автоматизированных систем. Научную деятельность 
Р. А. Нейдорф начал в НПИ (ЮРГТУ, г. Новочеркасск) как специалист по 
автоматизации химико-технологических процессов и продолжил в РГАСХМЕе и 
Донском государственном техническом университете (ДГТУ). 

Главной областью научных интересов Р. А. Нейдорфа была теория 


автоматического управления — в частности, разработка и развитие методологии 





структурного синтеза нелинейных систем автоматического управления. Это 
сложное направление на стыке интуиции, искусства и математики, требующее нестандартных идей и научной 
смелости в их отстаивании. Исследованиям в данной сфере посвящены кандидатская и докторская диссертации 
Р. А. Нейдорфа, монографии, публикации и доклады, а также диссертации многих его учеников. В широкий 
круг научных интересов Р. А. Нейдорфа входили также задачи и методы математического моделирования, 
оптимизации, обработки информации. Эти темы Р. А. Нейдорф разрабатывал со своими аспирантами. Итоги 
совместных исследований зафиксированы в публикациях и учебных пособиях. 

Результаты теоретических и опытно-конструкторских работ Рудольфа Анатольевича апробировались 
на крупных российских и зарубежных научных конференциях, внедрялись в производство. Один из примеров 
— метод экспериментального построения математических моделей существенно нелинейных систем, 
предложенный Р. А. Нейдорфом и группой молодых ученых. Решение нашло применение в выборе элементов 
структуры транспортной системы МААТ (МиБоду Адуапсе АпзМр Юг Тгапзрог®) (2012—2015 гг.). На нем 
базируется научный проект по теоретическому обоснованию метода экспериментального построения 
математических моделей и внедрению его в научную практику по гранту Российского фонда фундаментальных 
исследований. 

Р. А. Нейдорф внес существенный вклад в развитие инженерного образования в России, 
совершенствование учебного процесса, подготовку технических и научных кадров высшей квалификации. В 
вузах и на предприятиях Ростовской области работают два подготовленных Рудольфом Анатольевичем доктора 
наук и более 20 кандидатов наук в области математического моделирования, теории управления, методов 
исследования операций. 

Р. А. Нейдорф организовал в ДГТУ подготовку инженерных кадров по специальности «Управление и 
информатика в технических системах», научных кадров по специальностям «Математическое обеспечение и 
администрирование информационных систем» и «Программная инженерия». 

В 2002-2005 и 2006-2013 гг. Р. А. Нейдорф руководил кафедрами «Информационные и управляющие 
системы» в РГАСХМЕ и «Программное обеспечение вычислительной техники и автоматизированных систем» в 
ДГТУ. Он внес большой вклад в научную работу высших учебных заведений, позже объединившихся в 
региональный опорный вуз — ДГТУ, активизировал научно-исследовательскую работу в области 
информационных проблем управления. Р. А. Нейдорф, в частности, разрабатывал методы оптимизации, законы 
управления различными техническими системами, а также поисковые эвристические алгоритмы решения этих 
задач и установил приоритет исследовательских, опытно-конструкторских работ в данной области. 

В эти же годы он налаживал тесные связи с кафедрами «Управление в технических системах» 
Кисловодского гуманитарно-технического института и «Автоматизированные информационные системы» 


Ухтинского государственного технического университета, организовывал стажировки сотрудников. Последнее 
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время он в качестве главного научного сотрудника проводил большую работу в Центре научных компетенций, 
создавая методики подготовки кадров высшей квалификации. 

С 1996 года Р. А. Нейдорф возглавлял специализированный диссертационный совет Д.212.058.02. В 
2010 году он организовал ставший ежегодным семинар-конференцию «Системный анализ, управление и 
обработка информации». Рудольф Анатольевич был членом двух редакционных коллегий: журнала «Вестник 
Донского государственного технического университета» и научного сборника трудов организованного им 
семинара-конференции. Р. А. Нейдорф — автор и соавтор более 500 научных публикаций, в том числе 5 
монографий (3 зарубежные), 12 учебников и учебных пособий, более 70 патентов и авторских свидетельств. 

Не многим известно, что Рудольф Анатольевич был поэтом и редактором ежемесячной газеты 


«Атаманские вести», автором 158 поэтических произведений. 


Ищу мечту, а подбираю... прах, 
И потому нежней и осторожней 
Держу синицу — ту, что есть в руках. 
Так, говорят, вернее и надежней. 


И эту мудрость горькую мою 
Не упрекнут — вздохнут и оправдают. 
Петь в клетке тошно даже воробью, 
Подрезанные крылья не взлетают... 


И все же до конца безвестных дней, 
В туманной мгле призывный клин теряя, 
Благословляю в небе журавлей, 
Недостижимость их благословляю! 


Рудольф Нейдорф 
1984 год 
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Поперечные колебания круглого биморфа с пьезоэлектрическим 
и пьезомагнитным слоями 


А. Н. Соловьев, До Тхань Бинь, О. Н. Лесняк 
ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 


Введение. Статья посвящена исследованию поперечных осесимметричных колебаний биморфа с двумя 
пьезоактивными слоями: пьезоэлектрическим и пьезомагнитным. Этот элемент может найти применение в 
устройстве накопления энергии, которое находится в переменном магнитном поле. Целью работы являлось 
исследование зависимости частот резонанса, антирезонанса и коэффициента электромеханической связи от 
геометрических параметров элемента. 

Материалы и методы. Математической моделью работы пьезоэлемента является краевая задача линейной 
магнитоэлектроупругости. Элемент состоит из трех слоев: двух пьезоактивных (РИТ-4 и СоЕе›О4) и среднего 
пассивного слоя, выполненного из стали. В качестве метода решения краевой задачи используется метод 
конечных элементов, реализованный в пакете АМЗУ5. 

Результаты исследования. Разработана конечно-элементная модель пьезоэлемента в пакете АМЗУ$. Решены 
задачи определения собственных частот резонанса и антирезонанса. Построены графические зависимости этих 
частот и коэффициента электромеханической связи в зависимости от геометрических параметров устройства — 
толщины и радиуса пьезоактивных слоев. 

Обсуждение и заключения. Полученные результаты могут быть использованы в процессе конструирования 
рабочего элемента устройства накопления энергии за счет действия переменного магнитного поля. 
Построенные зависимости собственных частот резонанса и антирезонанса от геометрических параметров 
пьезоэлемента позволяют выбрать размеры пьезоактивных слоев для заданной рабочей частоты при 
наибольшем коэффициенте электромеханической связи. 


Ключевые слова: устройство накопления энергии, пьезоэлектрики, пьезомагнетики, метод конечных элементов, 
собственные частоты колебаний. 
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Введение. В создании сенсорных и измерительных систем, современной малогабаритной бытовой 
техники, телефонов сотовой связи и беспроводных сенсорных систем для наблюдения и диагностики 
технического состояния различных объектов не требуются мощные источники энергии, однако мобильность и 
энергонезависимость вышеперечисленных устройств являются обязательными. Устройства накопления энергии 
с пьезоактивными элементами, прямо преобразующими энергию механических колебаний в электрическую 
энергию, широко применяются для питания такого рода аппаратов. В [1-3] исследованы устройства накопления 
энергии с использованием пьезоэлектрических генераторов при воздействии механических нагрузок. 

Если система находится в переменном магнитном поле, созданном постоянными магнитами, 
установленными на вращающихся деталях машины, то пьезомагнитный слой деформируется вместе с 
пьезоэлектрическим элементом. За счет этого происходит образование электрического тока. Сио-Глап® Уа и др. 
в [7,8] обсуждали теоретические модели многослойных магнитоэлектрических композитов для 
магнитоэлектрического отклика на резонансных частотах, соответствующих колебаниям изгибной и 
экстенсиональной моды. Теоретическая модель магнитных энергетических собирателей, использующих 
функционально-градиентный композитный кантилевер, проанализирована для улучшения собирающей 
способности и регулировки резонансной частоты в [9]. В [10] У.Е. 7Бапга и др. изучали бифуркации, 
периодическую и хаотическую динамику четырехгранной композитной многослойной пьезоэлектрической 
прямоугольной пластины с простыми опорами. В [11] изучались изгибные и свободные колебания 
магнитоэлектроупругой пластины с поверхностными эффектами. 

В настоящей работе рассматриваются осесимметричные колебания устройства. Изучено влияние 
радиуса и толщины пьезоэлектрической и пьезомагнитной пластин на частотные характеристики устройства 
(собственные частоты резонанса и антирезонанса) и на эффективность преобразования энергии колебаний в 
электрическую энергию, которая характеризуется коэффициентом электромеханической связи. В качестве 
метода решения задачи выбран метод конечных элементов, реализованный в пакете АМЗУ$. 
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Рис. 1. Схема половины осевого сечения пьезомагнитоэлектрического генератора 
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Рассматриваемое устройство накопления энергии представляет собой осесимметричную конструкцию, 
которая состоит из металлического диска (подложка) с наклеенными на него двумя пьезоактивными слоями 
(рис. 1). Верхний слой —Й пьезоэлектрический, нижний — пьезомагнитный. Плоские поверхности 
пьезоэлектрического слоя покрыты электродами, которые подключены к внешней электрической цепи или 
один из электродов свободен, а на другом задан нулевой электрический потенциал. На пьезомагнитный слой 
действует переменное магнитное поле по гармоническому закону, внешний радиус подложки шарнирно 
закреплен. Математической моделью поперечных установившихся колебаний описанной конструкции является 
краевая задача линейной теории пьезомагнитоэлектроупругости [4]. В качестве инструмента решения выбран 
пакет АМЗУ5, в котором реализована модель пьезоэлектричества. 

Краевая задача для пьезомагнитоэлектрического тела состоит из системы уравнений и граничных 
условий [4]: 


У-с+рГ=рй, У.Р=а, У-В=0 (1) 
о=с:в-е'-Е-й'.Н 
Р=е:=+к:Е+а.Н (2) 
В=Й:=+а' -Е+и:Н 
1 
в=5(Ун +(Уи)'), Е=-Уф, в=-У$, (3) 
где о и = — тензоры механических напряжений и деформации; О и Е — векторы электрической индукции и 


напряженности электрического поля; Ви Н — векторы магнитной индукции и напряженности магнитного поля; 
р — плотность материала; с — тензор упругих модулей; е — тензор пьезоэлектрических модулей; й — тензор 


пьезомагнитных модулей; к — тензор диэлектрических проницаемостей; & — тензор магнитоэлектрических 
модулей; и — тензор магнитных проницаемостей; Г — вектор плотности массовых сил; с, — объемная 
плотность электрических зарядов; и — вектор перемещений; ф и ф — электрический и магнитный потенциалы. 
Граничные условия задаются для механических, электрических и магнитных компонент полей. 
Механические граничные условия. Пусть поверхность 5 состоит из двух частей Г, и Г,, так что 
5=Г, Г, , причём Г, ОГ, =. 
и=И0 на Г,, п:б=р на Г.. (4) 
Электрические граничные условия. Пусть поверхность 5 состоит из двух частей Г, и Г,, так что 


5=Г, 9 Г, , причём Г. ОГ, =©. 


ф=Ф(х,0) на Г,, пр =-в, на Г., (5) 
где о, — плотность поверхностных зарядов. Кроме того, в случае если электроды подключены к внешней 
цепи, необходимо добавить два условия: 

9]; =, |". Раб =Г, (6) 
5Е 
где 5, — площадь электрода; у — неизвестный потенциал, который находится из второго условия; / — 


электрический ток. 

Магнитные граничные условия. Пусть поверхность 5 состоит из двух частей Г; и Г., так что 
5=Г,; 9 Г, причём ГОГ, =©. 

ф= $(х,1) на Г;, п.В=ов, на Г, (7) 

где о, — плотность поверхностных свободных токов вдоль границы. 

Для упругого слоя используются первые уравнения в системе (1)-(3), неизвестными являются 
компоненты вектора смещений и, для пьезоэлектрического к ним добавляются вторые уравнения, неизвестны 
и и электрический потенциал ф, а для пьезомагнитного к первым уравнениям добавляются третьи и к 


неизвестным смещениям — магнитный потенциал ф. При этом соотношения (2) преобразуются обнулением 


соответствующих констант. 
Конечно-элементное моделирование. Компьютерная модель устройства построена в конечно- 
элементном пакете АМ5У$. Металлическая подложка (сталь) имеет толщину В и радиус г. Пъезоактивные слои 
состоят из одной пьезокерамической и одной пьезомагнитной пластин, которые поляризованы по толщине, 
имеют толщину В, и радиус г, (рис. 1). 
В конечно-элементной модели устройства для металлического и пьезоэлектрического слоев 
используются в качестве конечных элементов РГАМЕ42? и РГАМЕ1З3 соответственно. В этой работе 
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пьезомагнитный слой моделируется конечным элементом РГАМЕ1З3, в котором пьезоэлектрические свойства 
материала заменены на пьезомагнитные. Это возможно сделать по двум причинам: пьезослои между собой не 
контактируют, качественно уравнения для электрического и магнитного потенциалов совпадают. 

Свойства пьезоактивных слоев, используемые в расчетах, представлены в таблицах 1-2: пьезокерамика 
— РИТ-4, пьезомагнитный материал — СоЕе›О. [5, 6], клеевые слои не учитываются. 








Таблина 1 
Материальные свойства пьезокерамики РЯТ-4 
и, Гпа 2 Гпа ( ,Гпа Та | С. Тиа _ бзз, в а ие 
Кл/м” Кл/м” Кл/м” 
139 77,8 74,3 115 25,6 —5,2 15,1 12,7 730 635 
































Здесь 5) = 8,85х107'? Ф/м, плотность пьезокерамики РИТ-4 р = 7500 кг/м‘. 











Таблица 2 
Материальные свойства пьезомагнитного элемента СоЕе›Ол 
в ей ть т. тм О, Оз, О, Аз Азз› 
Гпа | Гпа | Гпа Гпа Гпа Н/Ам Н/Ам Н/Ам Нс’/Кл? Нс”/Кл” 
286 173 170 269,5 45,3 580,3 699,7 550 5,9х10 1,57х107" 
































Плотность пьезомагнитного элемента СоЕе›О.д р = 5290 кг/М?. 


Упругие свойства изотропного материала подложки характеризуются модулем Юнга ЕЁ и 
коэффициентом Пуссона У, в расчетах использовалиь Е=200ГПа, у=0,29, плотность р = 7860 кг/м? 


(сталь). 
Для достижения высокой точности в расчетах размер конечного элемента металлической слоя 
установлен значением не выше 1/5 его толщины, размер конечного элемента пьезослоев автоматически 


установлен. Сетка конечных элементов устройства накопления энергии представлена на рис. 2. 

















Рис. 2. Сетка конечных элементов устройства накопления энергии. 
Поверхности 1, 4 — внешние поверхности пьезослоев; поверхности 
2, 3 — контактные поверхности между пьезослоями и металлическим слоем 


При анализе собственных колебаний пьезоэлектромагнитного устройства накопления энергии 
предполагается выполнение следующих механических и электромагнитных граничных условий. На левом 
конце заданы условия симметрии, правый конец шарнирно закреплен (рис. 1). Для вычисления частоты 
резонанса устройства электрические потенциалы заданы на поверхностях 1, 2; на поверхности 3 задан 
магнитный потенциал; плотность магнитного потока задана на поверхности 4. В случае вычисления частоты 
антирезонанса устройства на поверхностях 1—4 заданы такие же граничные условия, как для вычисления 
частоты резонанса. Однако электрический потенциал на поверхности 1 неизвестен и находится из условия (6). 

Численные результаты. Рассматриваются собственные колебания пьезоэлемента, радиус которого 
составляет г,-9,8 мм, толщины пьезослоев В,-0,5 мм, радиус подложки г=10 мм, толщина подложки В=0,1 мм. В 
таблице 3 представлены первые три собственные частоты антирезонанса. 
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Таблица 3 
Собственные частоты антирезонанса 
Частота Распределение вертикального смещения на деформированном состоянии и 
№ | антирезонанса, распределение электрического и магнитного потенциала 
кГц на недеформированном образце 

р 

1 11,675 =— 
мх 

а Ш = в 
3 115,572 

у 

а 











В таблице 3 представлены распределения вертикального смещения, электрического и магнитного 
потенциала, при этом последние представлены в единой шкале, поэтому отображение магнитного потенциала 
не является наглядным. Для построения приближенной теории расчета колебаний рассмотренного биморфа, 
проведен анализ напряженно-деформированного состояния, электрического и магнитного полей, который 
показывает, что для низших изгибных мод могут быть принятые гипотезы об их распределении, 
соответствующие изгибу пластин. 

В представленных ниже результатах численных расчетов исследуется зависимость собственных частот 
резонанса и антирезонанса, коэффициента электромеханической связи от геометрических параметров. 

Когда число чередующихся пьезоэлектрических и пьезомагнитных слоев становится достаточно 
большим, возможно использование подхода на основе эффективных свойств пьезомагнитоэлектрического 
композита [5, 6]. В этом случае используются все уравнения задачи (1-3). 


Значение толщины пьезослоев В, изменяется в пределах 0,3-0,7 мм, а радиуса г, — в пределах 6,8-9,8 мм. 
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Рис. 3. Собственные частоты: а — частота резонанса, 
б — частота антирезонанса 


На рис. 3 представлены зависимости собственной частоты от значения толщины В, и радиуса г» 
пьезослоев. Из рис. 3 видно, что значения собственных частот возрастают с увеличением радиуса. 

На рис. 4 представлена зависимость коэффициента электромеханической связи от значения толщины В, 
и радиуса т, пьезослоев. Из рис. 4 видно, что значение собственной частоты возрастает с ростом радиуса 
пьезослоев г,, но убывает с ростом толщины пьезослоев В. 





Рис. 4. Коэффициент электромеханической связи 


Заключение. В работе рассматривается осесимметричная конечно-элементная модель устройства 
накопления энергии на основе круглых пластин в пакете АМЗУ5З. Активными элементами устройства 
накопления энергии являются пьезоэлектрическая и пьезомагнитная пластины, закреплённые на металлической 
пластине. Исследовано влияние геометрических характеристик пьезослоев в условиях определенных 
ограничений на размеры устройства, на коэффициент электромеханической связи, который показывает 
эффективность устройства накопления энергии. Результаты расчетов, приведенные в работе, позволяют 
выбрать рациональные размеры пьезоэлементов, работающих на определенной частоте и имеющие 
наибольшую эффективность. Данный анализ показывает, что максимальное значение коэффициента 
электромеханической связи достигается, когда толщина и радиус пьезослоев соответственно принимают 
наибольшее и наименьшее значения в рассматриваемых пределах. 
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Исследование распространения трещины в поверхностном белом 
слое рельсовой стали 


А. Ю. Перелыгина, В. Ю. Конюхов, А. Е. Балановский 
ФГБОУ ВО «Иркутский национальный исследовательский технический университет» (г. Иркутск, Российская Федерация) 


Введение. Статья посвящена оценке растрескивания белых слоев, формирующихся на поверхности рельса в 
процессе эксплуатации. Трещины в белом слое рельсовой стали обнаруживаются уже после 1 тыс. циклов 
испытаний. Это объясняется растягивающими и сдвиговыми напряжениями на поверхности пятна контакта 
колеса и рельса. В работе представлены результаты исследований морфологических характеристик белого слоя 
на поверхности рельса. 

Материалы и методы. Объект исследования (поверхность рельса после эксплуатации) рассматривался под 
микроскопом. Затем была разработана двумерная модель конечных элементов плоской деформации, чтобы 
моделировать динамические характеристики растрескивания белого слоя. Предложены математические модели, 
описывающие распространение трещины. Для этого применили критерий механики упругого пластического 
разрушения, метод /-интеграла. В программе 5У5ТЕГО выполнено численное моделирование образования 
белого слоя и распределения остаточных напряжений. 

Результаты исследования. Представлены оптические снимки микроструктуры поперечного сечения белого 
слоя на поверхности рельса после эксплуатации. На задней кромке белого слоя исследованных образцов 
зафиксированы два разных типа трещин. В программе 5У5ТЕГО визуализированы фрагменты моделирования 
механизма образования белого слоя с распределением остаточных напряжения, сжатия и растяжения. 
Результаты расчетов показывают, что значения ./-интеграла для всех трех трещин немного уменьшаются, если 
длина трещины достигает 10—50 мкм. 

Обсуждение и заключения. Полученные результаты применимы для оценки износостойкости рельсовых сталей 
и прогнозирования направления роста трещины. Сравнения /-интегральных максимумов показали, что при 
одинаковых условиях нагрузки трещина № 1, скорее всего, будет расти быстрее, чем трещины №2 из. С 
увеличением длины трещины максимумы ./-интеграла всех трех трещин уменьшались. 


Ключевые слова: рельсовая сталь, белый слой, трещина, упругое пластическое разрушение, метод /-интеграла 
распределением остаточных напряжений. 
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Введение. Известно, что детали машин и механизмов, различные функциональные конструкции часто 
испытывают катастрофическое хрупкое разрушение. При этом в зависимости от условий эксплуатации металл 
может разрушаться пластично или хрупко [1-3]. Хрупкое разрушение возникает из-за роста трещин, которые 
внезапно становятся неустойчивыми и распространяются в материале со скоростью звука. Трещины в металле 
могут иметь технологическое происхождение или зарождаться и расти в процессе эксплуатации. Белый 
нетравящийся слой (англ. уййепопесйтауег, ТЕГ — это явление, которое возникает на поверхности 
эксплуатируемого рельса под действием колес. И’ЕЁ образуется из-за сильной многоцикловой пластической 
деформации. Установлено [4—6], что от усталости при контактном качении на поверхности рельса трещины 
связываются с белым слоем. Слой назван так из-за его устойчивости к травлению кислотами во время 
металлографической подготовки и белого «безликого» внешнего вида (так он выглядит под микроскопом). 
Белый слой обычно встречается в пятнах контакта на поверхности рельса. Его глубина — -0,10-100 мкм после 
нескольких месяцев эксплуатации рельса [4]. Важная особенность белого слоя — твердость. Она достигает 
1300 НУ [5], но обычно находится в диапазоне 700-1200 НУ [6]. Такое состояние может привести к 
образованию хрупких трещин в белом слое и последующему распространению усталостных трещин. 
Авторы [4], исследовав рельсы после эксплуатации, обнаружили, что трещины присутствовали в белом слое, но 
не проникали за границу раздела с основным металлом рельса [4, 5]. Трещины развивались на поверхности и 
распространялись вдоль границы раздела между ЕЁ и перлитной структурой материала, вызывая сильный 
износ. Механизм образования белого слоя и его микроструктура изучались многими авторами [6-12]. 
Микроструктура слоя идентифицируется как перечисленные ниже смеси: 

— мартенсит с остаточным аустенитом и чистый мартенсит, насыщенный азотом и углеродом; 

— феррит с мартенситом; 

— остаточный аустенит и мартенсит с отпущенным мартенситом; 

— феррит с цементитом, карбид и нанокристаллическая фаза а-Ре. 

Обнаружено [6, 7, 9], что белый слой имеет мартенситную микроструктуру с высокой плотностью 
дислокаций. Авторы [8—10] предположили, что ЕЁ состоит из нанокристаллической фазы 9-Ре с размером 
зерен от 15 до 500 нм. Согласно [7], ЕЁ состоит из сильно деформированного перлита, нанокристаллического 
мартенсита, аустенита и цементита. В [8] представлены результаты рентгеновских измерений остаточных 
сжимающих напряжений в ЕЁ и выявлены значительные напряжения сжатия рельса (--600 МПа) в обоих 
направлениях: продольном и поперечном. Анализ работ [4-19] показывает, что трещины в белом слое 
рельсовой стали обнаруживаются после всех испытаний, даже тех, которые выполнены всего за | тыс. циклов. 
Причина тому — смешанные растягивающие и сдвиговые напряжения на поверхности пятна контакта колеса и 


рельса. После появления трещины в белом слое она быстро растет из-за хрупкости структуры мартенсита, пока 
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не достигнет границы раздела с основным металлом рельса, имеющим структуру перлита — здесь иная 
микроструктура отклоняет трещину и замедляет ее распространение. Это отклонение вызвано ориентацией 
цементитных пластин в перлитной структуре [10-19], которые параллельны поверхности качения. 
Согласно [12], снижение усталостного сопротивления рельса связано с давлением контакта и коэффициентом 
проскальзывания. Авторы [13] установили, что зарождение трещины в белом слое вызвано силой тяги и 
спроектированной поперечной нагрузкой сдвига к поверхности рельсовой колеи. Идентифицированы два типа 
трещин: 

— ведущая трещина дефекта в пятне контакта (вызвана сдвигом); 

— задняя трещина (хрупкое разрушение в форме клина). 

Авторы [14] связали коротковолновое динамическое взаимодействие между колесом и рельсом с 
началом и ростом деформирования поверхности. В условиях контактного качения в рельсовой стали при 
усталости выявили три вида трещин [15]: 

1) перпендикулярные поверхности контакта, 

2) с углом наклона, 

3) параллельные поверхности контакта на различных глубинах. 

Для определения напряженного состояния вблизи вершины трещины в условиях контакта с качением 
авторы предложили численный метод и зафиксировали изменение формы трещины. В [16] применили критерии 
направления роста трещины в смешанном режиме для оценки поведения при растрескивании. Показано, что 
трещины росли в направлении плоскости максимального напряжения сдвига, а не перпендикулярно плоскости. 
Авторы [17] разработали двумерную модель конечных элементов для моделирования поведения трещин в 
контакте качения (колеса и рельса) с учетом четырех коротких трещин. Обнаружено, что напряжение 
сдвигового режима играет доминирующую роль в росте трещины, более длинная трещина растет вверх, и это 
приводит к отколу поверхностного слоя. С помощью двумерной вычислительной модели авторы [18] 
обнаружили, что скорость роста трещин на поверхности рельса увеличивается с увеличением длины трещины и 
начинает уменьшаться после определенной глубины. В исследовании динамического взаимодействия колеса с 
рельсом [19] использована трехмерная модель конечных элементов и зафиксировано различное поведение при 
растрескивании из-за статических и динамических решений. В [20] учтены эффекты остаточного напряжения 
при изучении пути и скорости роста трещины. Во многих работах, в том числе [21, 22], показано влияние 
жидкости на дальнейшее развитие трещин. Предполагается, что жидкость на поверхностях трещин 
существенно влияет на их раскрытие (коэффициент интенсивности напряжений режима 1). Если же трещина 
развивается в сухих условиях, то доминирует режим сдвига [22]. До настоящего времени линейная механика 
упругого разрушения широко применялась при численном моделировании контакта колес с рельсами [23, 24]. 
Такой подход оценивает сопротивление нестабильному распространению трещины, или трещиностойкость 
металлов, при статическом нагружении по одному или нескольким критериям разрушения: 

— силовому (критический коэффициент интенсивности напряжений К 1 с); 

— деформационному (критическое раскрытие в вершине трещины); 

—щ энергетическому (критическое значение „Л-интеграла, работы пластической деформации и 
разрушения). 

На практике материалы вблизи области контакта могут быть пластически деформированы из-за 
высокой нагрузки на оси [3-6]. Чтобы получить более точное решение, следует использовать расчеты с 
упругопластической механикой разрушения или критериями, основанными на энергии. Кроме того, трещины 
могут формироваться в областях, прилегающих к белому слою [7-12], и механизмы их развития до сих пор 
полностью не изучены. Опубликованы работы по образованию белого слоя и трещин [14, 15], однако детально 
не исследовалось поведение распространения трещин в площади белого слоя. Более того, применительно к 
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взаимодействия. В данной работе представлены результаты исследования морфологических характеристик 
белого слоя на поверхности рельса. 

Материалы и методы. Сначала объект изучения рассматривался под микроскопом. Затем была 
разработана двумерная модель конечных элементов плоской деформации, чтобы моделировать динамические 
характеристики растрескивания белого слоя. Согласно выдвинутой гипотезе, область контакта на поверхности 
рельса пластически деформировалась при взаимодействии колеса и рельса, поэтому для оценки поведения 
распространения трещины введен /-интегральный критерий, основанный на энергии. 

Образцы, содержащие белые слои на поверхности головки рельса, вырезаны из рельса, в свое время 
эксплуатировавшегося на Восточно-Сибирской железной дороге. Рельс изготовлен из перлитной стали по 
ГОСТ51685-2013. Это сталь марки 76ХФ: С 0,78 %, 510,54 % и Ми 0,8 %, Ст0,40 %, У 0,035 %. Определены 
механические свойства рельса: 

— предел текучести — 944 кН, 

— временное сопротивление — 1287 кН, 

— относительное удлинение — 11,5 %, 

— относительное сужение — 31 %, 

— ударная вязкость 18 Дж/см”. 

Температура испытаний +20 °С. Макротвердость по Бринеллю на поверхности катания в местах без 
белых слоев составила 435, 485, 445, 465 и 494. Это не соответствует требованиям ГОСТ 51685—2013 (должна 
быть 352-405 по Бринеллю). Упрочнение поверхности головки рельсов в процессе эксплуатации обусловлено 
совокупным действием ряда физических механизмов [4-16]. Для расчета принят модуль Юнга Ё = 206 ГПа, 
модуль сдвига — 80 ГПа, плотность — 7850 кг/м? [3-10]. Стандартная металлографическая процедура 
выполнялась после секционирования образцов в направлении, параллельном направлению движения. После 
резки образец залили в смолу, прошлифовали с помощью шкурки карбида кремния (зернистость 1200-3000) и 
отполировали микросиликой 0,5 мкм. Поперечное сечение наблюдали в оптический микроскоп после травления 
2-процентным спиртовым раствором азотной кислоты в течение 7 секунд. Механизм образования белого слоя 
исследовали при помощи программы Б5УУТЕГО, которая использует аналитическую модель объемного 
тепловыделения двойного эллипсоида Голдака. Предписаны три типа теплового источника: 20 Гауссов, 
двойной эллипсоид и ЗО конический Гауссов. После того как определены геометрические параметры модели 
нагрева и задана величина максимального объемного тепловыделения, выполняется приближенный расчет при 
постоянных теплофизических характеристиках. Все необходимые для моделирования свойства материала 
задаются в виде кусочно-линейных функций. В программе 5УУТЕГ одновременно с тепловой задачей 
решается металлургическая. В последнем случае расчеты основаны на модели Леблонда, которая описывает 
металлургическое превращение одной фазы в другую. Решение задач механики сводится к определению 
термических деформаций. Для этого в программе нужно задать: модуль упругости первого рода, коэффициент 
Пуассона, коэффициент теплового расширения, а также кривые упрочнения для металлургических фаз. При 
растрескивании в белом слое в условиях подвижного контакта наблюдения выполняли с использованием 
программного обеспечения (ПО) АМ№М5У5/ [5-РУМА, которое предназначено для моделирования методом 
конечных элементов. При исследованиях динамики железнодорожных экипажей широко применяется 
известный алгоритм названного ПО — ЕА5Т$/М с эллиптической областью контакта [25]. По сравнению с 
геометрическими размерами колеса и рельса область контакта довольно мала (площадь не более 15 мм’ [25]. 
Следовательно, модель контакта колеса с рельсом можно упростить до условия деформации 20 
плоскости [26, 27]. На продольном сечении были построены белые слои в форме дуги длиной 2 мм и толщиной 
0,2 мм (параметры получены в результате анализа исследований И’ЕЁ на поверхности рельса [4-—15]). С целью 
снижения вычислительных затрат и обеспечения точности в модели применили метод адаптивного построения 


сетки. В районе белого слоя использовались очень мелкие сетки рис. 1. 
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Рис. 1. Конечноэлементное моделирование (КЭМ) белого слоя с тремя трещинами (а); и преобразование интеграла контура 
элемента в интеграл контура поверхности: произвольный контур вокруг вершины трещины (6); различные контуры 
элементов, выбранные для представления независимости от пути вычисления ./-интеграла (в); принципиальная схема 20 
метода виртуального закрытия трещины для расчета коэффициента интенсивности напряжений (2). 


Установленный минимальный размер элемента — 0,02х0,02 мм”. Элементы на концах и гранях трещин 
имели более высокую плотность сетки с установленным минимальным размером элемента 0,01х0,01 мм”. Это 
позволяет точно захватывать градиенты напряжений. За пределами зоны растрескивания применены более 
крупные сетки с минимальным размером элемента 0,1х0,1 мм” в контактных поверхностях колеса и рельса. Они 
постепенно увеличиваются по направлению к дальним участкам поля. Сначала рельс объединили в 776 919 
элементов с 789 084 узлами, которые изменились во время моделирования из-за применения метода адаптивной 
сетки. Такой размер модели является приемлемым для получения точных расчетов контакта с качением [19, 20, 
28]. Постоянное вертикальное усилие нагружения в 13000 Н было приложено к колесу, что соответствовало 
эквивалентному максимальному давлению по Герцу -1,2 ГПа [25]. Колесу задали скорость вращения ® 
43,5 рад/с, что равняется скорости перемещения 72 км/ч с коэффициентом трения 0,3. Исследовано, как на 
свойства трещин в белом слое влияют их длина, изменения угла наклона, давление нагрузки и коэффициент 
трения. 

Рост трещин, разрушающих поверхность, проанализирован в предположении линейной механики 
упругого разрушения [1-3, 17-19]. Для более точного решения использован критерий механики упругого 
пластического разрушения, метод /Линтеграла, который применим как для линейных, так и для упругих 
пластических решений [29]. Он был введен для изучения поведения распространения трещины. Для численных 
расчетов /-интеграл можно получить из решения в дальней зоне [29]. В рамках работы с конечными 
элементами интеграл контура элемента должен быть преобразован в интеграл контура поверхности, как 
показано на рис. 16, где Г — это контур кривой, который определяет границу ./-интеграла, и он направлен 
против часовой стрелки от нижней грани трещины к верхней. /-интеграл может быть дополнительно оценен 
как: 


1 = (мау-Т, <), 
где » — плотность энергии деформации; Т; — вектор тяги; и; — компонент вектора смещения; 4 — 
приращение длины по контуру Г [29]. 
Оценка ./-интеграла реализована в [/5-РгеРо5! как инструмент постобработки. Это приложение также 
подходит для 2) моделирования плоской деформации [30]. Для расчета „/-интеграла вокруг вершины трещины 
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были выбраны различные контуры элементов (рис. 16). Важно отметить, что условия разгрузки для 
независимого от пути /-интеграла не предлагаются в соответствии с механикой разрушения [31]. Принцип 
расширения Л-интегральной методологии за пределы допустимости линейной механики упругого разрушения 
заключался в идеализации упругопластической деформации как нелинейной упругой деформации. Известна 
идентичность характеристик нагрузки — деформации для упругопластических и нелинейно-упругих 
материалов [30, 31]. Упругопластический материал следует линейному пути разгрузки с наклоном, равным 
модулю Юнга. Нелинейный упругий материал разгружается по тому же пути, что и при загрузке. Таким 
образом, анализ предполагает, что нелинейно-упругое поведение справедливо для упругопластического 
материала. Разгрузка происходит вблизи вершины трещины после того, как колесо прошло по рельсу. 
Рассчитанные значения ./-интеграла могут быть действительными только до достижения максимума. Поэтому в 
представленной работе обсуждаются только максимальные показатели „-интеграла. Коэффициенты 
интенсивности напряжений в рамках конечноэлементной схемы вычислялись 20 методом виртуального 
закрытия трещины [29—32]. Силы у и х, необходимые для соединения узлов си 4 (рис. 1г), обозначаются Р’. и То 
соответственно. 

Следует отметить, что теория механики линейного упругого разрушения основана на предположении 
об отсутствии пластической деформации вокруг трещин. В то же время, согласно работам [1-3], использование 
теории механики линейного упругого разрушения корректно, т. к. нелинейная деформация материала 
ограничена небольшой областью, окружающей кончик трещины. Некорректны рассуждения о том, что в связи с 
малой областью контакта можно рассматривать плоскую задачу. При плоской деформации эта область 
полосовая и бесконечная. А в рассматриваемом случае речь идет о существенно трехмерной задаче. В [32] 
описывается сложность геометрии рельса и ограниченность области, в которой происходят основные 
изменения напряженно-деформированного состояния. Авторы названной работы предлагают рассматривать 
упрощенную осесимметричную задачу для многослойного покрытия в форме кусочно-неоднородного слоя. 
В [33] исследуется напряженно деформированное состояние многослойного покрытия в окрестности области 
контакта колеса и рельса. При этом показано, что жесткость и толщина верхнего слоя покрытия существенно 
влияют на эквивалентные и контактные напряжения. При фиксированных механических параметрах покрытия 
с увеличением толщины его верхнего слоя уменьшаются значения максимальных эквивалентных и контактных 
напряжений. Авторы [34] наносили тонкие многослойные покрытия на поверхность железнодорожного рельса 
в окрестности бокового контакта рельса и колеса. Напряженно-деформированное состояние этих покрытий 
исследовалось при различных значениях их геометрических и механических параметров в процессе 
поворотных движений. 

Результаты исследования. На рис. 2 представлены оптические снимки микроструктуры поперечного 
сечения белого слоя на поверхности рельса. 
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Рис. 2. Особенности трещин в белом слое на поверхности рельса, обнаруженные с помощью оптической микроскопии: 
общий вид белого и деформированного слоя при 200-кратном увеличении (а); 1 — трещина переднего края, 
2 — трещина среднего положения при 500-кратном увеличении (6) 


На рис. 2а показан сегмент белого слоя на поверхности. Его глубина 10 мкм. Значение микротвердости 
в белом слое НИ! = 989. Это соответствует особенностям ЕЕ, представленным в работах [4—15]. Толщина 
ТЕТ, может достигать 200—300 мкм в зависимости от условий эксплуатации [4—6]. На изображении видна резкая 
граница между белым слоем и деформированной перлитной структурой основного металла рельса. Глубина 
зоны пластической деформации между ТЕГ и основным материалом рельса составляет -70 мкм. По 
оценкам [5-10], на эксплуатируемых перлитных рельсах в пятне контакта «колесо — рельс» пластическая 
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деформация имеет высокую концентрацию, создавая зону тонкого сдвига 502 180 мкм. В то же время 
испытания на твердость показывают, что из-за объемного контакта глубина пластической деформации может 
доходить до 1-10 мм [4-9, 27—30]. Это приводит к важной эволюции микроструктуры в более узком масштабе, 
а именно: к изгибу или разрыву цементитных пластин, уменьшению межпластиночного расстояния. Меняются 
и механические свойства, т. к. увеличивается прочность на растяжение и уменьшается вязкость разрушения для 
трещин, параллельных выровненным пластинкам цементита. На переднем крае ИЕГ расположены две 
различные трещины (рис. 26). Они могут возникать и распространяться на границе раздела между ТЕГ и 
деформированной перлитной микроструктурой. Детальное исследование других образцов показало, что 
наклонная трещина появилась на передней кромке ЕЕ, следуя направлению течения материала рельса. Затем 
трещина идет вниз, в материал рельса вдоль границы слоя. На рис. 26 показана трещина, распространяющаяся 
вертикально в среднем положении в пределах ТЕГ. На исследованных образцах на задней кромке ЕЁ 
наблюдаются два разных типа трещин. Они получены из разных участков ТЕГ, но расположены на одном и 
том же рельсовом образце и очень близко друг к другу. Трещина первого типа формируется и распространяется 
вдоль границы между белым слоем и зоной пластической деформации (рис. 26). Трещина второго типа 
пересекает границу раздела с деформированной микроструктурой, но демонстрирует тенденцию искривления 
вдоль линии по направлению деформации сдвига перлита. 

Такое поведение трещин описано в работах [4—7]. Действительно, исследованные нами трещины 
расположены в одном и том же белом слое, следовательно, они подвергались одинаковым нагрузкам. Однако 
скорости их роста и развития различны. Это указывает на изменение состояния напряжения в ведущем, среднем 
и заднем положении в белом слое. Значит, должен быть другим механизм последующего распространения 
трещины (т. е. позади ее). Как показано в работах [35, 36], пластически деформированная структура перлитной 
области непосредственно под белым слоем может играть существенную роль в воздействии на распространение 
трещин. Трещины, которые образуются из-за усталости при контакте качения, можно разделить на две 
категории по месту их образования: под поверхностью и на поверхности. 

Обычно под поверхностью трещины возникают из-за сильной вертикальной нагрузки в сочетании с 
дефектами материала. На поверхности большинство трещин образуются из-за взаимодействия колеса и рельса, 
а также переноса большой нагрузки в небольшую область контакта. Площадь контакта эллиптическая. Она 
относительно мала и при этом поддерживает всю нагрузку колеса. Трещины, которые образуются в результате 
усталости при контакте качения и являются результатом интенсивного напряжения сдвига в области контакта 
колеса и рельса, будут расти, когда эти напряжения превышают величину допустимого растяжения рельсовой 
стали. Также возможно, что трещины продвинутся к вершине рельсов. Согласно [11, 12, 23], скорость нагрева 
поверхностного слоя рельса, возникающая при прохождении поезда, может превышать 10°°С/с и достигать 
температуры, при которой образуется аустенит. После прохождения поезда при быстром охлаждении аустенит 
превращается в мартенсит. Время контакта колеса с рельсом чрезвычайно коротко (миллисекунды). При этом 
сталь из эвтектоидного состава нагревается от комнатной температуры до 727°С со скоростью 10°°С/с. Таким 
образом, практически невозможно измерить фактические изменения температуры в течение рассматриваемого 
периода. С другой стороны, быстрый нагрев и закалку можно моделировать и контролировать в программе 
УУ5ТЕГО. Для этого использовали скорости нагрева 20-1000 °С/с, которые получали за счет параметра 
скорости перемещения теплового источника по поверхности рельса в диапазоне скоростей 5-—100 мм/с. 

На рис. 3 показан фрагмент моделирования механизма образования белого слоя в программе 5У5ТЕГО. 
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Рис. 3. Результаты моделирования формирования белого слоя на поверхности головки рельса в программе ЗУЗ\У/ЕГО: доля 


мартенсита в белом слое, % (а); распределение остаточных напряжений в белом слое, МПа (6) 
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По данным табл. 1 можно судить о том, как уровень остаточных напряжений зависит от параметров 
дорожки и скорости упрочнения. 
Таблица 1 
Величина остаточных напряжений (сжатия и растяжения) по ширине и глубине дорожки при скорости 
упрочнения 15 мм 





Расстояние по центру оси 
дорожки (глубина), мм 
Остаточные напряжения, МПа —43 —87 —15 323 445 347 396 274 176 
Расстояние по центру от оси 
дорожки (ширина), мм 
Остаточные напряжения, МПа —43 —92 | -140 | -189 | -116 30 200 249 128 


0 0,65 | 0,77 | 0,86 | 0,96 | 1,22 | 1,54 2,5 3,0 








0 0,14 | 1,65 | 1,92 | 2,19 | 2,59 | 3,16 | 3,91 5,2 









































Белый слой, смоделированный в программе, состоит из мартенсита. Морфология слоя аналогична 
наблюдаемой в исследуемом рельсе (рис. 2). Твердость моделируемого ИЕГ, — 670-810 НТ. Это немного ниже, 
чем в рельсах (ТЕГ - 725-1050 НГ). Однако текущие моделирования представляют только один цикл фазового 
превращения. Моделировались от одного до пяти повторных циклов термообработки на поверхности головки 
рельса. Установлено, что твердость увеличивается до 700-850 НУ. Очевидно, мартенсит с такой же твердостью, 
что и в белом слое рельса, может быть получен после нескольких проходов колес. При моделировании 
формирования белого слоя выяснилось, что в этом процессе важную роль играют изменение температуры во 
время контакта колеса с рельсом и распределение остаточных напряжений. 

На рис. 4. показаны /-интегральные максимумы на кончиках смоделированной трещины. 
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Рис. 4. Максимальные значения интеграла на кончиках трещин в передней, средней и задней кромках ЕЁ 
(длина трещины для расчета выбирается как 0,01 мм; 0,03 мм; 0,05 мм; 0,1 мм и 0,2 мм отдельно) 


Следует отметить, что длина трещины измеряется от поверхности рельса на глубину ИЕГ. На 
представленном варианте моделирования трещины № 1,3 возникают под углом а ир (83,13°), максимальное 
нагрузочное давление — 1,2 ГПа. При этом все три трещины имеют длину 10 мкм. На вершине трещины № 1 
значение ./-интеграла достигает максимума 0,321 Н / мм. Максимальные значения ./-интеграла для трещин № 2 
и №3 намного меньше — 0,18 Н /мм и 0,088 Н / мм соответственно. Анализ рис. 4 показывает, что трещина 
№ 1 на переднем крае 7ЁГ, скорее всего, будет расти по сравнению с трещинами в среднем положении и на 
задней кромке. Этот факт согласуется с экспериментами на усталость при качении [10-18]. Показано, что 
граница (начало) белого слоя на поверхности рельса и основного металла имеет самую низкую 
сопротивляемость образованию дефекта контактно-усталостного характера, в то время как усталостная 
долговечность в конце ЕЁ примерно в 3 раза выше, чем в центре и начальной точке. Из рис. 4 видно, что 
максимумы -интеграла вокруг вершин трещин показывают тенденцию к уменьшению с увеличением длины 
или глубины. Это явление более значительно для трещины № 1, чем для № 2 и № 3. По мере увеличения длины 
трещины максимальное значение Линтеграла для всех трех трещин уменьшается соответственно с 
0,3214 Н/мм до 0,222 Н /мм; с 0,18 Н /мм до 0,14 Н/мм ис 0,088 Н / мм до 0,062 Н / мм. Это согласуется с 
результатами [19], поскольку на более глубокие поверхности трещины меныше влияют нормальные и 
тангенциальные контактные нагрузки. Однако из рис. 4 видно, что для трещины № 2, когда ее длина составляет 
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200 мкм и соответствует толщине ЕЕ, максимум /-интеграла демонстрирует противоположную тенденцию, о 
которой никогда не сообщалось. Максимум ./-интеграла в этом случае наибольший по сравнению с другими 
короткими трещинами № 2. Можно предположить, что трещина № 2 будет быстро распространяться к границе 
МЕТ, и подложки. Это хорошо видно на рис. 26 (стрелка 2). Существенная разница объясняется материальными 
неоднородностями между идеально упругой ТЕЁ и упругопластичной рельсовой матрицей. Результаты 
расчетов показывают, что для всех трех трещин значения Линтеграла немного уменьшаются, когда длина 
трещины увеличивается от 10 до 50 мкм. Если углубление трещины №2 составляет от50 до 100 мкм, 
изменение Л-интеграла относительно невелико. Напротив, отмечается существенный рост максимального 
значения /-интеграла, если трещина № 2 распространяется через весь слой на границу между ЕЁ и матрицей 
рельса. 

Максимальное значение ./-интеграла резко возрастает с усилением нагрузочного давления от 0,8 
до 1,5 ГПа: 

— с 0,27 Н /мм до 0,53 Н /мм для трещины № 1, 

— с 0,17 Н /мм до 0,47 Н /мм для трещины № 2, 

— с 0,04 Н / мм до 0,1 Н /мм для трещины № 3. 

Это говорит о том, что трещины с большей вероятностью распространяются при высоких нагрузках. 
Максимальное основное напряжение вокруг вершины трещины можно использовать как индикатор для оценки 
свойств разрушения и пути трещины [19-26]. Например, если длина трещины №1 составляет 10 мкм, 
максимальное основное напряжение на вершине трещины — 1589 МПа. Для трещин №2,3 с одинаковой 
длиной получаются разные концентрации напряжений вокруг вершин трещины, и максимальное основное 
напряжение составляет около 1506 МПа и 1261 МПа соответственно. Увеличение длины с 10 до 100 мкм 
приводит к снижению максимального основного напряжения с 1506 МПа до 1206 МПа. Когда трещина 
распространяется вниз к границе белого слоя, можно наблюдать резкий рост концентрации напряжений. Это 
хорошо видно из моделирования распределения остаточных напряжений по глубине слоя (рис. 3). Итак, 
согласно „/-интегральным расчетам, при одинаковых условиях нагрузки трещины в передней кромке белого 
слоя распространяются с большей вероятностью, чем в средней позиции и задней кромке. Это согласуется с 
ростом трещин. При этом из-за неоднородности материала возникает сложное поле напряжений, когда средняя 
трещина распространяется на границе раздела белого слоя и основного металла. Это приводит к высокой 
концентрации напряжений и тем самым ускоряет распространение трещины. Следует отметить, что в данной 
работе не учитывались анизотропные характеристики материала, вызванные сильной пластической 
деформацией в переходной зоне. В [4-10] показано, что такого рода деформация сдвига приводит к удлинению 
и выравниванию цементитных колоний поверхностного слоя головки рельса (рис. 2а). Выявлено, что полосы 
сдвига и области с высокой плотностью дефектов встречаются в переходной зоне белого слоя и некоторые 
цементитовые пластинки в этой области разрушены [4-8]. Хотя механика линейного упругого разрушения не 
отражает детальную физику ситуации [1-3], она способствует пониманию направлений наблюдаемого роста 
трещины. Для получения более развернутой информации по динамике трещинообразования в белом слое на 
поверхности рельса необходимо учитывать коэффициенты интенсивности напряжений. В то же время трещины 
контакта качения испытывают смешанную нагрузку, соответствующую последовательности растягивающих 
напряжений с последующим циклом сдвига [4—20]. Следовательно, коэффициенты интенсивности напряжений 
в направлении максимального тангенциального напряжения также могут применяться для прогнозирования 
направления роста трещины. 

Обсуждение и заключение. Установлено, что трещины в переднем крае белого слоя (трещина № 1) 
распространяются от края слоя вдоль границы раздела между белым слоем и основным металлом рельса. 
Трещины №2 (в середине слоя) растут вертикально, прежде чем достигнуть пластически деформированной 
матрицы. Трещины № 3 задней кромки демонстрируют две тенденции в развитии. В первом случае трещина 
идет по границе раздела, во втором пересекает материал рельса и постепенно выравнивается по направлению 
деформации. Сравнения интегральных максимумов показали, что при одинаковых условиях нагрузки 
трещина № | с наибольшей вероятностью будет иметь высокую скорость роста по сравнению с трещинами 
№2, 3. С увеличением длины трещины максимумы ./-интеграла всех трех трещин уменьшались. Установлено 
большое влияние давления нагрузки на развитие всех рассмотренных типов трещин. 
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МАШИНОСТРОЕНИЕ И МАШИНОВЕДЕНИЕ 
МАСШ\Е ВОПЛИ С АХО МАСЧИТ\Е 5СТЕМСЕ 








УДК 519.6 Врз://401.0г2/10.23947/1992-5980-2020-20-2-137-142 
Устойчивость нелинейно-упругой пластины при боковом сжатии [=] =] 

Г. И. Волокитин 

ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) [=] 


Введение. Если круглая пластина нагружена по боковой поверхности, возможны потеря ее устойчивости и 
выпучивание. Данную проблему можно решить с учетом бифуркации. В этом случае пластина рассматривается 
как нелинейно-упругое тело. В задачах устойчивости нелинейной упругости важен выбор зависимости между 
напряжениями и деформациями. В ранних работах, посвященных этой проблеме, рассматривались простые 
законы состояния (конститутивные уравнения). В качестве примера можно привести материал «гармонического 
типа» Сенсенига. 

Материалы и методы. Для круглой пластины из материалов Мурнагана и Блейтца и Ко получены уравнения 
нейтрального равновесия. В предположении однородной начальной деформации рассмотрена задача 
устойчивости этой пластины. Строгие трехмерные уравнения нейтрального равновесия позволяют исследовать 
смежные формы равновесия с учетом физической и геометрической нелинейности. Вывод этих уравнений 
основан на теории наложения малой деформации на конечную. 

Результаты исследования. Продвижение в решении соответствующего векового уравнения (с нелинейным 
вхождением параметра) для практически важных законов упругости Мурнагана и Блейтца и Ко возможно лишь 
с использованием численных методов. Разработанный метод расчета бифуркационных значений параметров 
нагрузки дает возможность проанализировать влияние нелинейности. 

Обсуждение и заключение. Исследовано влияние физической и геометрической нелинейности на величину 
верхнего критического значения параметра начальной деформации. Полученные результаты могут быть 
использованы при оценке достоверности модулей упругости третьего порядка для разных физических 
материалов. Данные об этих модулях пока малочисленны. Численные исследования показали, что к 
приведенным в некоторых источниках константам следует относиться с осторожностью. Обсуждается также 
использование модулей упругости в законе состояния Блейтца и Ко. 


Ключевые слова: конечная деформация, напряжения, выпучивание, пластина, устойчивость. 
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За ШиИу оЁГа попПтеаг @а$Яс р! \е ипдег ]жега] соттргез$10п 
С. Г. Удокт 


Доп Зе Тесиса! ОшуегзИу (Козюу-оп-Роп, Кл531ап Еедегайоп) 


тодисйоп. Г.о$$ ОЕ За бИИу апа БисКПпе оРа гоипа рае тау Бе обзегуе4 1 Фе р|=е 15 1оа4е4 оп ®е 1айега] зиГЁасе. 
'ТБе зоГайоп ю 1$ ргоМет 15 Базе оп а ЫЁигсайоп арргоасй. ш 1$ сазе, а рае 15 сопз14еге4 аз а попПпеаг е|а5с 
Боду. И 15 паронапЕ {0 своозе фе геайоп Бебмееп ге55е5 ап аеЮюгтаНопз ш зизатаб у ргоепл$ оР попПпеаг 
е]азнсиу. Зпаре 1а\уз оЁ зе (сопзвийуе едиайоп5) \еге сопз1Аеге ш еаЙу \огК$ деусе4 1ю 15 ргоет, юг 
ехатр[е, таепа| оЁ Ве “Баглотс фуре” (Зепзеп12). 

Маепаб апа Мейо%. Едлайопз оР пеига1 едигиит ог гоипа р]аез таде оР МигпазВап апа В1ай-Ко табепа1$ 
аге оМалте4. Аззитте а ипопи пива! дегоппайоп оп фе ре, фе за БШиу ргоет 1$ сопз1егед. Зилс{ гее- 
Аитепз1опа| пешга] еда Гат едиайоп$ рго\1е ехрогие геже4 огз оР едаЬгат 1аКше ш® ассоипЕ рБуз1са1 
ап4 сеотей1с попПпеагИу. ОепуаНоп оЁ Фезе едиаНоп$ 1$ Базе оп Ше аррИсаНоп оЁ Ше еогу оЁ зарегроз оп оЁ 
эта дегогтайоп оп фе Йпа] опе. 
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ппроцапе 1а\уз оЁ @азйсйу оЁ МигпазВап апа В!а-Ко 1$ роз51е изте фе питетса! тефо4$ оту. ТБе 4еуе!ореа 
тео4 ог са!сшайпе ЫРагсайоп уаез о оа@ те рагатеегз таКез { роз е ю апа[уте 1е еЁЙес{ оР попПпеагиу. 
Р5сиз$1от апа Сопсия1опя. ТВе шНчепсе оЁ рБуз1са| ап4 зеотейлс попПпеагиу оп Фе иррег см@са| уае оЁ Фе 
ша] 4еГоглайоп рагатеег 1$ ехроге4. ТВе гезиНз оатед сап Бе изе4 ипаег Пе аззеззтлеп{ оЁтепабИиу оЁ е!а5ис 
{фиа-отдег тодиЙ Юг уапоч$ рБуз1са|! таёепа15. аа оп Фезе то@иН 15 $ИП зсагсе. ТБе питепса| гезеагсв Ваз 
зво\уп фай 1е сопзвапз отуеп ш зоте зоигсез зношШ@ Бе шеа{е4 у саиноп. Те изе оРе]азНсИу пода ш Ше 1а\\у оЁ 
бабе оЁ В1а-Ко 1$ а|50 415си$$еа. 


Кеуптога$: Ипие деогтаНоп, з#геззез, БискПпо, р!айе, аби. 


Еог сйайоп: С. 1. Уо]окит. ЗабИйу оРа попПпеаг е!азйс ре ипаег ]а{ега| сотргез$1оп. УезиК о ОЗТЦ, 2020, 
у01. 20, по. 2, рр. 137-142. Б&рз://401.ог/10.23947/1992-5980-2020-20-2-137-142 


Введение. В настоящее время изучение новых сравнительно простых и адекватных законов состояния 
различных материалов, требующих учета нелинейных деформаций, представляет собой активно развивающееся 
направление механики сплошной среды [1-10]. В рамках теории наложения малой деформации на конечную 
выведены трехмерные уравнения нейтрального равновесия для материалов Мурнагана и Блейтца и Ко. На 
основе этих уравнений рассмотрен пример торцевого выпучивания равномерно сжатой круглой 
пластины [11, 12]. 

Материалы и методы 

Уравнения нейтрального равновесия. Пусть г, 9, 2 — цилиндрические координаты точки в 
недеформированном состоянии. Предполагаем, что начальная деформация тела определяется радиус-вектором 
В: 

В =А(ке, + 421, (1) 


где К(г) = аг ‚ параметры а, 4 — константы, е,, е» в — базисные векторы. 


Для координат этой точки в начально-деформированном состоянии имеем: 
В=аг, Ф=о Г=ф. 
Следовательно, градиент деформации, мера деформации Фингера и ее главные инварианты определены 
соотношениями: 
рт — ., — хр2 
УК=УЕ =а(е,е, +ерер) +4, Е =УК`, 


п =2а’+4*, Г =а+2а’4`, [=а\а?. (2) 
Здесь и далее У, У — набла-оператор в метрике недеформированного и начально-деформированного 


состояния: У =УЮ`".У. 


Будем использовать уравнения нейтрального равновесия, предложенные А. И. Лурье [13]: У@ =0, где 
тензор @ — линейный дифференциальный оператор над вектором добавочных перемещений У’. Выражения 
компонент этого тензора, конкретизированные с учетом законов состояния Мурнагана и Блейтца и Ко, 
получены в работах [14, 15]. В представлениях тензора © компоненты — некоторые функции, определяемые в 
результате решения краевой задачи начальной деформации. 

Результаты исследования. Для рассмотрения изгибной формы бифуркации равновесия пластины 
аналогично работам [11, 12] принимаем следующий вектор добавочного перемещения 

УИ! = и(г,2)е, + и(г,2)4 . (3) 
Учитывая (1), (2) и (3), тензор Лурье имеет вид: 
9 = Де, е, + Ве е, +Сы, + Се, + Нед. 
Здесь 
да а И вы Е Е Пао 
г д» 02 г д» 02 г д 





о 
02 д’ 02 Ге) 


Отметим, что структура этих операторов типична в задачах устойчивости цилиндрических нелинейно- 
упругих тел. Опуская выражения остальных коэффициентов, приведем, например, формулы для А. В случае 





закона состояния Мурнагана [13, 14] А выражается соотношением 


“(2 Зу, +4, 





А = 


У 
Е : четы 
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Гоомит С. 1. МЧаБИйу ора потПтеаг ейа5йс р/ще ипаег [ег сотргезоп 





Здесь и далее ^. и и — модули упругости Ляме, у, у,, у, — константы упругости третьего порядка. 
Для материала Блейтца и Ко [15]: 
А и 2и(1 =. В) ра + В аа , 
ат В 


где В — уточняющий модуль упругости. 
Уравнение нейтрального равновесия равносильно системе дифференциальных уравнений относительно 


компонент вектора У’: 
2 1 2 2 
а [2 м [до вобу-о 








РИ а дубё а ' 02? 4 
2: 2 
ыы), о и 
02\ 0" г д гб’ а ` 02 


Предполагая, что по боковой поверхности пластина нагружена равномерным давлением, 
дифференциальные уравнения (4) дополним условиями равновесия на границе: 





и(",2)]_ =0, ©9,.|_, =0. 
На торцах, т. е. при 
р 
а 9.=0, 9. =0. (5) 
Применим подстановку 
и=Х, (2). Л (КГ), 6) 
и= 2, (2).Ло (А„г). 
Здесь п=1,2,..., /ь, /, — функции Бесселя нулевого и первого порядка, а числа К,’, — нули функции Л,(х). 
Примем х = Аг и применим равенства 
[9 4?.Л (К АЛ, (Е 
м ло, АО -юло, О ло --ьло, 
а’ Ч’ Ч’ 


а также известные тождества для функций Бесселя 
1 1 1 
п ' п" ' 
лы, ИЕЛыЕЛ, 
х х х 


Выполнив перечисленные выше действия, получим, что переменные в дифференциальных уравнениях 


(4) и граничных условиях (5) разделяются. Приходим к следующей краевой задаче для обыкновенных 
дифференциальных уравнений: 





А 
Хх, г г р К, С —_ 7, = 0, 
Н, С 
(7) 
[е; 6,+С, 


ЕЯ, +, 


п 





Х! =0. 
3 3 
СХ, (+#/2)+С,2, (+1/2)=0, (8) 
К,6,2, (+1/2)-С,Х! (+1/2) =0. 
Два однородных линейных обыкновенных дифференциальных уравнения (7) и четыре краевых 
условия (8) приводят к задаче на собственные значения с нелинейным вхождением параметра. В данной задаче 
таким параметром является величина (1 _ а) . Систему (7) приведем к нормальному виду: 


У =У,», 


А О. +С 
= уча, — у,, 
Н, \@р 





1 
Уз = У, 
С С, +С 
' 2 1 2 2 
= К, Уз К, У. 
С, С, 
Т Т 
Здесь приняты обозначения: аргумент 1=ЕФ, (х; у; У.) = (х,; ЖИ 2) . Пусть 
> Т 
фундаментальная система решений — это следующие четыре вектора у, =(у,;; у: У; У,). Тогда общее 
Т 4 
решение системы имеет вид (Х,; ЛХ’; 7,; 2) и = >, у, где & — произвольные постоянные. Зададим, 
= 


например, при 2=-й/2 начальные данные столбцами единичной матрицы четвертого порядка. Затем решим 
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численно задачу Коши с этими начальными условиями. В результате в правой части получим значения 
базисных функций, то есть векторы У;(+ й/2). Используя граничные условия (8), приходим к однородной 


системе линейных алгебраических уравнений относительно &, &,, 6, &: 
4 
а, =0, где 1=1, 2, 3, 4. 
= 


Здесь коэффициенты — элементы матрицы А: 
а, =ЮСу, (-1/2) + СУ. (-1/2) ‚а, = СУ (1/2) + Су (1/2) › 


а; =К, С, Уз; (-#/2)- СУ, (-1/2) ‚ 44, =К, СУ, (й/2) — б1у,; (1/2) у 
Однородная система имеет нетривиальное решение, если выполнено условие 
4е/А=0. (9) 
Выражение определителя включает параметры нагрузки а, 4, нули функции Бесселя, а также модули 


упругости Л, и, у, у,, ъ, (для материала Мурнагана) либо /., и, В (для материала Блейтца и Ко). 


12 


Параметры а, @, задающие начальную деформацию, взаимосвязаны. Осевая сила, действующая на 
К, 

площадку поперечного сечения, определяется соотношением [13]: О=2л | о,КаК, где о, — физическая 
0 


компонента тензора напряжений. Поскольку начальная деформация предполагается в виде (1), то тензор 
напряжения Коши Т и мера деформации Фингера Г соосны. Причем тензор напряжений — постоянный: 


2 г: ы . г 
Т = а" (а (ее, +ее,)+4°ъ)-с(а“(е,е, +е,е,) +4“) +сС'(е,е, +ее, + ьь,)) : 


Торцы пластины свободны от нагрузки, поэтому о, =0. Отсюда следует условие, связывающее 


коэффициенты закона Фингера 
а 
с® с +-——=0. (10) 
а? 


Принимая во внимание закон состояния, можно записать конкретное выражение условия (10), которое 
устанавливает связь между аи 4. Так, для материала Мурнагана получаем условие: 


ре +[еи- Зи —4у, + зи - 








3 15 
3).—2 + “. + —- +, + (2^-Зу, —4у, _ + (у +у, — а“ =0. 
В частности, если у, =у, =у, =0 ‚то имеем 
4? Ри. 
А+2и 
При рассмотрении закона состояния Блейтца и Ко 4 выражается формулой 
ей 


Таким образом, задавая значения модулей упругости, нули функции Бесселя, из (9) находим 
бифуркационные значения параметра начальной деформации а,. 


Отметим, величина Л =|-а, для сравнительно тонких дисков — малый параметр. Так, в классической 
теории выпучивания пластин критическое значение величины а, определяется формулой 
2 
в (3,8317) [ й 
12+) (и, 


где у — коэффициент Пуассона, а 3,8317 — первый корень функции Бесселя .Л(х). 





* 


При решении начальной краевой задачи предполагается возможность применения метода возмущений 
Синьорини. При этом считаем, что коэффициенты оператора ® зависят от малого параметра Л степенным 
образом. Это означает, что условия на границе для компонент добавочного перемещения задаются Л матрицей 
А= А(Л). В данной задаче имеет физический смысл частичная проблема собственных значений — 
определение низших собственных чисел. Старшие степени Л, которые содержат элементы матрицы А, 
незначительно сказываются на значениях наименьших корней векового уравнения (9). Если при решении 
начальной задачи ограничиться линейной теорией, то А(Л) — регулярный двучлен [16]. 

Численные эксперименты при исследовании устойчивости нелинейно-упругих тел не слишком 
большой относительной толщины подтверждают этот вывод. Следовательно, характеристическое уравнение с 
нелинейным вхождением параметра (9) можно заменить характеристическим уравнением линейного оператора. 


Волокитин Г. И. Устойчивость нелинейно-упругой пластины при боковом сжатии 


Гоомит С. 1. МЧаБИйу орга поттеаг ейа5йс р/ще ипаег [ег сотргезз?оп 





Можно применить итерационные процессы, сходящиеся к одному собственному числу, где в качестве нулевого 
приближения выбиралось число, близкое к значению в теории пластин. 

Обсуждение и заключение. С помощью полученных выше уравнений выполнен численный анализ 
влияния физической и геометрической нелинейности на величину верхнего критического значения параметра 
а,. Расчеты реализованы в среде МаЙаЬ для различных вариантов задания модулей упругости, относительной 


толщины пластины и числа волнообразования п. Обнаружено, что во всех случаях потери устойчивости диска 


минимальному критическому значению параметра а отвечает первая форма выпучивания, которая 


соответствует корню функции Бесселя 3,8317 [17]. 
В табл. 1 приведены критические значения параметра а, для пластин относительной толщины 0,05; 0,1; 
0,15; 0,2; 0,25 и 0,3. 




















Таблица 1 
Критические значения параметра а, для пластин различной относительной толщины 
в у Относительная толщина пластины 
и 0,05 0,1 0,15 0,2 0,25 0,3 
| 0,9922 0,9869 0,9800 0,9660 0,9530 0,9330 
0,9972 0,9867 0,9735 0,9549 0,9300 0,8914 
2 0,9977 0,9908 0,9801 0,9664 0,9504 0,9332 
0,9967 0,9869 0,9690 0,9410 — —- 
3 0,9976 0,9907 0,9800 0,9662 0,9502 0,9330 
0,9965 0,9856 0,9665 0,9372 — — 
я 0,9977 0,9907 0,9798 0,9659 0,9497 0,9324 
0,9881 — — 0,9965 0,9613 0,9861 
5 0,9914 0,9865 0,9783 0,9668 0,9520 0,9350 
0,9929 0,9729 0,9341 0,9114 —- — 
6 0,9978 0,9912 0,9800 0,9675 0,9520 0,9340 
0,9985 0,9941 0,9871 0,9775 0,9651 0,9539 
7 0,9985 0,9941 0,9859 0,9735 0,9560 0,9320 
0,9985 0,9941 0,9863 0,9749 0,9594 0,9396 





























Здесь цифрами обозначены перечисленные ниже материалы. 

Сталь Вех 535 (\.=1,09, и=0,818, у =-1,75 , у, =-2,40 ‚ у, =-1 69). 
Сталь 50ХГСМ2Ф (^. = 1,129, и = 0,803, у, =-2,8 , у, =-2,7, у, =-187). 
Сталь Неба 37 (^. =1,11, и=0,821, у, =-3,58, у, =-2,82, у,=-177). 
Сталь Неа АТУ (^=0,87, и=0,716, у, =0,34, у, =-5,52, у, =-Ь0). 
Бронза бериллиевая (^,= 1,042, и=0,49, у, =-4,0, у, =-17, у, =-0,6). 
Оргстекло (^=0,39, и=0,186, у, =-0,078, ъ, =-0,07, ъ, =0,047). 


О о 


В первых шести вариантах рассмотрен закон состояния Мурнагана. Последний вариант представляет 
результаты для материала Блейтца и Ко. Верхнее число в ячейке таблицы относится к случаю, в котором не 
учитывается физическая нелинейность, т.е. у =у, =у, =0 в законе Мурнагана; р =0 в законе Блейтца и Ко. 


Второе (нижнее) число учитывает физическую нелинейность. В законе Мурнагана были использованы данные 


дин 


> 12 
для модулей упругости третьего порядка из [14] в единицах 10 о 


В последней строке табл. 1 приведены результаты для закона Блейтца и Ко, который выбран в 
упрощенном варианте (уравнение Ноулса — Стернберга): коэффициент Пуассона принимался равным 0,25, а 
уточняющий модуль В = 0,45. Прочерк означает, что критическое значение а, не обнаружено. 

Анализ результатов позволяет сделать некоторые выводы. При малых относительных толщинах диска 
точная теория и линейная теория пластин дают одинаковые критические значения параметра а.. 
Геометрическая нелинейность заметно влияет при относительных толщинах больше 0,1. Еще сильнее 
сказывается физическая нелинейность. Однако следует с осторожностью выбирать значения упругих модулей 
третьего порядка. Например, в четвертом варианте (сталь Нес!а АТУ) и для тонких пластин не были 
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обнаружены надежные критические значения параметра а,, хотя для тех же значений модулей упругости Ляме 


наблюдаются приемлемые значения, свидетельствующие о потере устойчивости. Как и в задаче устойчивости 
нелинейно-упругой сферы из материала Блейтца и Ко [15], константа В слабо влияет на величину критического 


значения а,. 
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Разработка методики проектирования технологического процесса 
обработки шарико-стержневым упрочнителем с учетом 
формирования сжимающих остаточных напряжений 





М. А. Тамаркин, Э. Э. Тищенко, С. А. Новокрещенов, С. А. Морозов 
ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 


Введение. Представлены результаты исследований процесса обработки многоконтактным виброударным 
инструментом — шарико-стержневым упрочнителем. Описаны методы обработки поверхностным 
пластическим деформированием, их преимущества. Приведена схема инструмента. Для шарико-стержневого 
упрочнения выявлены его особенности, технологические преимущества, область применения. 

Материалы и методы. При проведении теоретических исследований процесса обработки установлены 
факторы, влияющие на качество поверхностного слоя обработанных деталей. Приведены зависимости для 
расчёта шероховатости поверхности, глубины упрочненного слоя, степени деформации при обработке шарико- 
стержневым упрочнителем. При исследовании формирования остаточных напряжений выявлена зависимость 
для расчета остаточных напряжений, образующихся в поверхностном слое обработанной детали. 

Результаты исследования. Представлены результаты экспериментальных исследований процесса обработки, 
необходимых для проверки адекватности приведённых теоретических моделей, а также методика их 
проведения. Приведены таблица и графики, наглядно подтверждающие хорошую сходимость теоретических и 
экспериментальных данных (разница не превышает 20%). Остаточные напряжения в поверхностном слое 
являются сжимающими, что позволяет прогнозировать высокие эксплуатационные свойства обработанных 
деталей. Величина остаточных напряжений на поверхности детали находится в пределах 130-200 МПа. 
Глубина залегания сжимающих остаточных напряжений находится в пределах 0,9-—1 мм. 

Произведены расчеты изменения усталостной характеристики — предельных напряжений цикла по глубине, 
которая оказывает влияние на величину предела выносливости. Установлено, что обработка деталей шарико- 
стержневым упрочнителем позволяет повысить предельное напряжение цикла при циклическом нагружении 
детали на 27-35%. 

Обсуждение и заключения. Предложенная методика проектирования технологического процесса обработки 
шарико-стержневым упрочнителем может быть использована при разработке технологии на 
машиностроительных предприятиях. В соответствии с рекомендациями задаются пределы необходимых 
параметров качества поверхностного слоя обрабатываемой детали, выбираются параметры шарико- 
стержневого упрочнителя, натяг и радиус заточки стержня. Производится расчет параметров качества 
поверхностного слоя. Корректировка выбранных режимов и повторный расчет параметров обработанной 
поверхности производятся до тех пор, пока все заданные характеристики не будут располагаться в 
необходимых пределах. 


Ключевые слова: обработка шарико-стержневым упрочнителем, шероховатость поверхности, глубина 
упрочненного слоя, степень деформации, остаточные напряжения. 
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Введение. Качество деталей зависит от множества различных параметров, так или иначе 
проявляющихся в процессе их изготовления и влияющих на их жизненный цикл. Особую роль в формировании 
качества поверхности занимают финишные методы обработки, использование которых позволяет улучшить 
основные эксплуатационные свойства деталей машин, такие как контактная жесткость, износостойкость, 
усталостная прочность и долговечность и т.п. 

Методы обработки поверхностным пластическим деформированием (ППД) позволяют решать задачи 
увеличения жизненного цикла изделия для деталей из различных материалов любой конфигурации. Кроме того, 
такие методы позволяют осуществлять обработку локально, что позволяет значительно снизить расходы на 
производство. Характерными поверхностями деталей, обрабатываемыми местным упрочнением ППД, 
являются: зоны концентрации напряжений (отверстия, шлицы, скосы, выборки, резьбы, галтели, пазы и т.д.); 
неупрочненные участки поверхности деталей, прошедших общую упрочняющую обработку в вибрационных, 
дробеструйных и прочих установках (под прижимами, в карманах, отверстиях и других труднодоступных для 
обрабатывающей среды зонах); места механической доработки деталей и др. 

Методы обработки ППД, в том числе и локальные, позволяют значительно снизить шероховатость 
поверхности обрабатываемой детали, повысить ее физико-механические свойства. Кроме того, методы ПИД 
позволяют снимать в детали остаточные напряжения, образованные различными дефектами, которые 
оказывают значительное влияние на физико-механические свойства материала, из которого деталь изготовлена. 
Известно, что растягивающие остаточные напряжения могут значительно снизить прочностные характеристики 
детали, и наоборот, сжимающие — их улучшить. 

Для осуществления процесса местной обработки ППД деталей простой и сложной конфигурации, в том 
числе имеющих небольшой перепад профиля по высоте, на кафедре «Технология машиностроения» Донского 
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государственного технического университета под руководством профессора Бабичева А. П. было изобретено 
специальное устройство — шарико-стержневой упрочнитель (ШСУ). ШСУ представляет собой 
многоконтактный виброударный инструмент, обработка которым построена на основе поверхностно- 
пластического деформирования. Данный метод обработки имеет ряд преимуществ, таких как возможность 
обработки концентраторов напряжений крупногабаритных изделий, изделий сложной конфигурации, 
маложестких изделий. Метод сочетает технологические возможности вибрационной отделочно-упрочняющей 
обработки и чеканки. 

В качестве силового привода используется пневмомолоток, на котором закреплен корпус упрочнителя. 
Боек силового привода | наносит удары с частотой порядка 42 Гц по нескольким слоям стальных шаров 4, 
которые передаются пакету круглых стержней 2, установленному в цанговом зажиме 6 (рис. 1). Наличие 
нескольких слоев шаров позволяет стержням, имеющим сферическую заточку рабочей поверхности, 
копировать фасонный профиль обрабатываемой детали 5. 
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Рис. 1. Схема многоконтактного виброударного инструмента ШСУ: 
1 — силовой привод, 2 — пакет круглых стержней, 3 — корпус упрочнителя, 
4 — стальные шары, 5 — обрабатываемая деталь, 6 — цанговый зажим, 7 — упругий элемент 


Основная часть. 

Рассматриваемый метод обработки исследован довольно большим количеством авторов [1-10]. При 
исследовании технологических возможностей процесса обработки ШСУ было установлено, что наибольшее 
влияние на качество поверхностного слоя обработанных деталей оказывает энергия удара привода, диаметр 
заточки стержней, число стержней в пакете, подача устройства вдоль обрабатываемой поверхности [1, 2]. В 
работе Щерба Л. М. [1] разработана методика проектирования технологического процесса обработки шарико- 
стержневым упрочнителем. Представлены зависимости для определения параметров качества поверхностного 
слоя: шероховатости обработанной поверхности, степени деформации и глубины упрочненного слоя. В работе 
Исаева А. Г. [2, 5] эти зависимости уточнены и описаны ниже. 

Шероховатость поверхности определена по формуле 


2 = 0,03 





(0 


где Е, — энергия удара индентора; № — число стержней в пакете; НВ — твердость обрабатываемого 





материала по Бринеллю; 1 — коэффициент полезного действия устройства, зависящий от ряда факторов (натяг 
при обработке, количество слоев шариков и т.п.), р — диаметр сферической заточки стержня упрочнителя 
(индентора). 

Значение степени деформации представлено как: 





(2) 





(3) 


В работе Копылова Ю.Р. [3] приведена зависимость определения остаточных напряжений для 
вибрационной отделочно-упрочняющей обработки. Представляется возможным применить аналогичную 
методику для определения остаточных напряжений при обработке ШСУ, введя коэффициент КА, ‚ учитывающий 





особенности формирования внутренних напряжений при обработке ШСУ. Коэффициент А, определен при 
проведении экспериментальных исследований. Зависимость после преобразований примет следующий вид: 
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Е, 1-У2 1-2 (3 
{е ЕАК. У - ы д т ш , (4) 
ра Е, в 
2(р+4) 
где К, = коэффициент, учитывающий особенности формирования внутренних напряжений при обработке 


ШСУ (определяется экспериментально), А 


‚о — коэффициент, учитывающий способность материала к 


упрочнению; В, = энергия удара индентора; Рида — диаметры заточки стержня и отпечатка на детали 


соответственно; У, и У, — коэффициенты Пуассона для детали и стержня; Ё, и Ё, — модули упругости 


первого рода для детали и стержня соответственно. 
Для проверки адекватности зависимости (4) 
формирования остаточных напряжений при ШСУ. Измерение остаточных напряжений осуществлялось по 


проведены экспериментальные исследования 
стандартной методике, используемой заводской лабораторией ПАО 
автоматизированного стенда контроля остаточных напряжений АСКОН-3-КИ на образцах алюминиевого 
сплава В95. Обработка прямоугольных образцов 200х100х20 мм осуществлялась на различных режимах с 


натягом 1,5 мм и 4,5 мм. Использовались пакеты стержней с радиусом сферической заточки 4 мм и 8 мм. 


«Роствертол» с применением 


Определение остаточных напряжений осуществлялось по методу Давиденкова путем стравливания 
поверхностных слоев с образцов-пластин, вырезанных из обрабатываемых прямоугольных образцов. 
Измеренные значения сравнивались с теоретическими, рассчитанными по формуле (4). Результаты расчетов 
представлены в таблице 1. 

Значение коэффициента К 


и 


определено по результатам предварительных экспериментальных 
исследований путем сравнения величины остаточных напряжений, рассчитанных По зависимости (4) и 


полученных экспериментально. Он изменяется в диапазоне 1,5-1,7. Для расчетов принято значение 1,6. 


Таблица 1 
Сравнение результатов теоретических и экспериментальных исследований 




















Разница теоретических 
Исходные данные 
Теоретические Экспериментальные данных от 
Радиус заточки т значения, МПа данные, МПа экспериментальных 
стержня данных, % 
4 4,5 141,6 168 15,7 
4 1,5 137,3 135 1,7 
8 4,5 1241 130 4,5 
8 1,5 118,74 126 5,6 

















На рис. 2,3 представлены графики зависимостей остаточных напряжений от натяга при обработке. 
Сплошной линией показана теоретическая кривая. Квадратами обозначены результаты экспериментальных 


исследований. 
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Рис. 2. Зависимость остаточных напряжений от натяга при обработке. 
Радиус сферической заточки стержня В=4 мм 
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Рис. 3. Зависимость остаточных напряжений от натяга при обработке. 
Радиус сферической заточки стержня В=8 мм 


Для полученных в результате исследований значений остаточных напряжений произведены расчеты 
изменения усталостной характеристики — предельных напряжений цикла по глубине по формуле, приведенной 
в работе [1]. Известно, что на величину предела выносливости оказывает влияние величина среднего 
напряжения цикла, которое в значительной степени зависит от уровня остаточных напряжений поверхностного 
слоя. 

Графики представлены на рис. 4-7. Сплошной линией показаны предельные напряжения цикла, 
пунктиром — распределение остаточных напряжений. 
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Рис. 4. Изменение величины предельных напряжений цикла Рис. 5. Изменение величины предельных напряжений 
по глубине детали после обработки ШСУ: К =4 мм, цикла по глубине детали после обработки ШСУ: В=4 мм, 
натяг — 4,5 мм натяг — 1,5 мм 
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Рис. 6. Изменение величины предельных напряжений цикла Рис. 7. Изменение величины предельных напряжений цикла 
по глубине детали после обработки ШСУ: К=8 мм, по глубине детали после обработки ШСУ: В=8 мм, 
натяг — 1,5 мм натяг — 4,5 мм 
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Заключение 

По результатам проведенных исследований можно сделать следующие выводы. 

1. Остаточные напряжения в поверхностном слое являются сжимающими, что позволяет прогнозировать 
высокие эксплуатационные свойства обработанных деталей. 

2. Величина остаточных напряжений на поверхности детали находится в пределах 130—200 МПа, что по 
данным ПАО «Роствертол» примерно соответствует величине остаточных напряжений после упрочняющей 
вибрационной обработки. 

3. Глубина залегания сжимающих остаточных напряжений находится в пределах 0,9-1 мм, что по данным 
ПАО «Роствертол» примерно в два раза больше, чем при упрочняющей вибрационной обработке. 

4. Обработка деталей ШСУ позволяет повысить предельное напряжение цикла при циклическом 
нагружении детали на 27-35%. 

5. Режимы обработки оказывают незначительное влияние на величину и глубину залегания остаточных 
напряжений. 

6. Установлено, что разница между результатами теоретических и экспериментальных исследований не 
превышает 20%, что свидетельствует об адекватности полученной теоретической зависимости для расчета 
остаточных напряжений при обработке шарико-стержневым упрочнителем. 

Исследования остаточных напряжений при ШСУ представлены в работах Шведовой А. С. [4, 6-9]. В работе 
[4] произведено компьютерное моделирование процесса с помощью пакета программ Апзуз. Полученные значения 
остаточных напряжений коррелируют с результатами вышеприведенных экспериментов. 

По результатам проведенных исследований предлагается следующая методика проектирования 
технологических процессов обработки шарико-стержневым упрочнителем. 

1. Задаются пределы необходимых параметров качества поверхностного слоя обрабатываемой детали. 

2. Производится выбор пневмолотка. Затем осуществляется выбор насадки с учетом имеющихся в 
распоряжении технолога. Для обработки малых участков рекомендуются насадки с малым числом стержней, для 
обработки участков большой площади — с болышим числом стержней. Число слоев шаров выбирается в 
зависимости от высоты перепадов либо радиуса кривизны обрабатываемой поверхности. 

3. Выбирается натяг при обработке и радиус заточки стержня. Для более твердых материалов подбираются 
большие значения натяга обработки и меньшие радиусы заточки. 

4. Назначаем время обработки 10-15 секунд на площадь пучка стержней. С учетом рекомендованного 
выше значения подачи выбирается число проходов инструмента по поверхности детали, при этом в большинстве 
случаев желательно использовать обработку в один проход. При наличии нескольких вариантов сочетаний режимов 
обработки, позволяющих получить заданные характеристики поверхностного слоя, выбирается тот, у которого 
общее время обработки конкретной детали меньше. 

5. Затем производится расчет среднего арифметического отклонения шероховатости обработанной 
поверхности, глубины упрочненного слоя, степени деформации, остаточных напряжений по формулам 1-4. 

6. По результатам расчетов производится корректировка выбранных режимов обработки. Затем вновь 
рассчитываются параметры обработанной поверхности и так до тех пор, пока все заданные характеристики не будут 
располагаться в необходимых пределах. 

Приведенная методика может быть использована при технологическом проектировании процесса 
обработки ШСУ. 

Вышеприведенные технологические рекомендации были использованы при внедрении технологического 
процесса обработки ШСУ на ПАО «Роствертол». 
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Моделирование интегрального показателя качества автомобиля [2] м] 
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ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) [=] 





Введение. Работа посвящена сравнению и оценке качества легковых автомобилей с использованием математи- 
ческого моделирования интегрального показателя качества. Цели работы: 

® разработка метода математического моделирования интегрального показателя качества; 

® составление поэтапной последовательности оценки качества; 

® анализ наиболее распространенных на отечественном рынке марок легковых автомобилей; 

® подведение итогов и обоснование полученного результата. 
Материалы и методы. Предложена последовательность моделирования интегрального показателя качества. 
Проведено численное моделирование на примере шести марок легковых автомобилей, наиболее распростра- 
ненных на отечественном рынке. 
Результаты исследования. Разработан метод моделирования интегрального показателя качества. Описана по- 
этапная последовательность проведения анализа. Предложена аддитивная форма определения интегрального 
показателя качества, позволяющая объединять в процессе сравнения единичные, не связанные между собой 
показатели качества в комплексный показатель. 
Обсуждение и заключения. Представленные предложения могут быть использованы при оценке и диагностике 
конкурентоспособности не только автомобилей, но и других товаров, обладающих большим набором независи- 
мых друг от друга показателей качества. 


Ключевые слова: моделирование, интегральный показатель, требования качества, автомобиль, технические 
характеристики, аддитивный метод, коэффициент весомости. 
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Введение. Для получения, систематизации, описания и применения наработанных знаний и инфор- 
мации в любой сфере деятельности наиболее широко применяется универсальный метод моделирования. Ма- 
тематическое моделирование систематически подкрепляется и актуализируется во многих областях наук. Ма- 
тематическое моделирование позволяет объединять наработанные знания и реальные процессы, в том числе 
мыслительные [1], объединять возможности воспроизведения свойств реального или создаваемого объекта, 
процесса или явления с применением другого объекта, процесса или явления. 

Качество является совокупностью потребительских свойств, не связанных между собой техническими 
характеристиками, которые должны удовлетворять текущим и перспективным потребностям [2]. При этом 
можно проследить тесную взаимосвязь качества с требованиями. Требования — это конкретные особенности 
и условия, соответствующие потребительским предпочтениям в течении всего срока эксплуатации. Требова- 
ния к товару закладываются еще на стадии проектирования [3, 4]. 


Однако, качество и требования имеют некоторое несоответствие [5, 6]. Особенно это проявляется в 
несоответствии заявленных технических качеств товара и требований, предъявляемых потребителями. Это объ- 
ясняется постоянным изменением требований потребителей, зависящих от технического прогресса, а также 


финансового и культурного состояния населения. 

В зависимости от количества характеризуемых свойств и от их влияния на качество изделия показа- 
тели качества подразделяются на единичные и комплексные. Единичный показатель позволяет охарактеризо- 
вать одно свойство или одну зависимость качества от технико-эксплуатационной характеристики, тогда как 
комплексный показатель объединяет несколько характеристик. Существует понятие прямых и обратных пока- 
зателей. Увеличение прямых показателей приводит к повышению качества, увеличение обратных показателей 
— кпонижению [7]. 

Проанализировав классификацию показателей качества продукции (рис. 1), можно делать вывод о 
том, что интегральный показатель является частным случаем комплексного показателя. Интегральный пока- 
затель качества характеризует отношение суммарного полезного эффекта при эксплуатации данного продукта 
к суммарным затратам на его приобретение, эксплуатацию или потребление [8]. 


Показатели качества продукции 
Единичные ЗЕЕ Комплексные 


Относитель- Опреде- 
| ее А 


Рис. 1. Классификация показателей качества продукции 










Метод моделирования интегрального показателя качества изделий. Определение качественного 
показателя на примере легкового автомобиля осуществляется в последовательности: 

® постановка цели; 

® подбор прототипов для сравнения и выбора оптимального варианта; 

® подбор параметров, характеризующих изделие с позиции потребителя; 

® определение усредненных значений показателей качества; 

® объединение усредненных данных в один показатель для каждой марки; 

® анализ полученных интегральных показателей всех прототипов; 

® принятие решений по управлению уровнем качества продукции. 

Наибольшие трудности представляют подбор характеристик и получение усредненных значений этих 
данных [9], так как нет известных методов формирования списка показателей качества и трансформации каж- 
дого показателя в числовую форму. Для определения интегрального показателя качества широко применяют 
его аддитивную форму, т. е. средневзвешенное суммирование, позволяющее объединять единичные, не свя- 
занные между собой показатели качества в комплексный показатель. Его применяют в случаях принятия ре- 
шения по экономически наиболее выгодному варианту, если на рынке представлено несколько прототипов. 
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То есть, данный метод является наиболее подходящим для моделирования интегрального показателя качества 


автомобиля. 
Формула для расчета комплексного интегрального показателя [10]: 
— уп у 
К: = 2-10: `ХЬ 
где а; — коэффициент весомости 1-го параметра; Х; — усредненный показатель качества по 1-му параметру; и 


— число параметров для сравнения. 
Для определения коэффициента весомости применяется экспертный метод в соответствии с выраже- 
нием: 


Алгоритм применяемых расчетов можно проследить на примере выбора оптимального варианта приоб- 
ретения легкового автомобиля (рис. 2). 
52 373 ил, 


| АДА бгапка ы | № 


39 149 шт. 
КА о 


24 687 шт, 
Нуллба1 $08 


22 395 ил; 
кокаови РоЮ 


14 312 шт, 
кода Вард 





Ве 
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Рис. 2. Популярные автомобили типа седан отечественного рынка 





Анализу подвергались легковые автомобили, наиболее продаваемые на отечественном рынке. Иерар- 
хия осуществлена по ценовой характеристике. 

Результаты исследования. Для чистоты эксперимента сравнивались шесть автомобилей с кузовом ти- 
па седан в одной ценовой категории. Показателей качества несколько десятков, в расчетах же рассчитан инте- 
гральный показатель для 6-ти групп основных (с точки зрения потребителя) показателей качества. Исходные 
данные для расчета приведены в табл. 1. На практике применения подобных оценок качества можно учитывать 
разное количество характеристик, даже с учетом того, что многие них имеют достаточно противоречивое влия- 
ние на качество и по своей физической сущности являются весьма сложными для их количественной оценки. 


























Таблица 1 
Исходные данные для расчета интегрального показателя качества автомобилей 
Модель 
о Гад Нуцп4а: | Код к: Уок 
ада уип4а1 оЧа та оКзуавеп 
6 Г. 
г Оташа Зо]а115 Кара Вю А Роо 
Цена, тыс. руб 580 750 680 640 600 690 
Пробег, тыс. км 50 75 57 26 53 50 
Расход топлива, л 7,2 8,3 8,5 6,4 8,7 8,0 
Мощность, кВт 72,07 90,50 80,91 73,50 78,00 80,91 
Клиренс, мм 170 160 136 160 ИТ 163 
Время разгона, с 11,5 11,0 10,0 13,1 12,8 11,7 





























Промежуточным этапом вычислений является получение усредненных значений 6-ти рассматриваемых 
показателей качества для каждой модели Коэффициент весомости определяется по методу экспертной оценки, 
где человеческий фактор (личностные предпочтения) также играет немаловажную роль при оценке и суще- 
ственно влияет на итоговое значение интегрального показателя качества сравниваемых седанов. Результаты 
расчета представлены в табл. 2. 


Таблица 2 
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Расчетные параметры качества автомобилей 























п Модель 
ры Ад ба Нуипаа1 ЗКо4а Ка Гаа УоЖК5\азеп 
о Зо1а11$ Вара Вю Уеза Роо 
Место по объему 1 й 6 3 2 5 
продаж 
К 
ОИ 0,11 0,15 0,25 0,15 0,14 0,20 
весомости 
р 0,32 0,32 0,71 0,48 0,33 0,49 
показатель качества 
Место по 5 5 1 3 4 2 
коэффициенту качества 
Конкурентоспособность 55:10 42-10 103 75:10 55:10 71-10 





























Таким образом, максимальный интегральный показатель качества принадлежит автомобилю ЗКо4а 
Кар! и составляет 0,71. Вместе с тем значения коэффициента качества остальных сравниваемых автомобилей, 
находящихся в одном ценовом диапазоне, отличаются незначительно. Значения коэффициента конкурентоспо- 
собности сравниваемых седанов получены отношением коэффициента качества к стоимости автомобиля. Этот 
показатель также максимален для Кода ВКар14, что объясняется тем, что покупатель данного автомобиля на 
одну денежную единицу приобретает больше качества, чем покупатели сравниваемых прототипов. 

Обсуждение и заключения. Фиксация информации об объекте исследования или проектирования для 
хранения и передачи в пространстве или во времени требует применения математического моделирования. Мо- 
делирование направлено на построение, совершенствование, изучение и применение моделей реально суще- 
ствующих или проектируемых объектов для последующей оценки их конкурентоспособности. Интегральный 
показатель качества является одним из унифицированных показателей, способных реализовать сравнение пока- 
зателей даже с противоречивым влиянием на качество. 
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Проектирование мобильных машин моделированием динамических 
нагрузок на узлах их приводов 





С. А. Партко, Л. М. Грошев, А. Н. Сиротенко 


ФГБОУ ВО «Донской государственный технический университет, Ростов-на-Дону, Российская Федерация» 


Введение. Проектирование мобильных машин невозможно без учета колебательных параметров их агрегатов, 
что требует разработки специализированных динамических моделей, учитывающих вероятностный характер 
этих параметров. Основная причина возникновения колебательных воздействий обусловлена неровностью 
профиля пути мобильной машины и непостоянством физико-механических характеристик почвы. Задачи, 
учитывающие эти особенности, с достаточной точностью решаются линейно, но в многомерных динамических 
системах такой подход недопустим ввиду наличия большого количества взаимных связей. Цель работы — 
сравнительный анализ эффективности существующих методик расчета статистических характеристик 
некоррелированных внешних воздействий применительно к мобильной машине, представленной в виде 
многомерной динамической системы с воздействиями, имеющими различные корреляционные связи. 
Материалы и методы. Внешние воздействия в многомерных динамических системах рассмотрены в 
матричной форме. При вычислении статистических характеристик учтены взаимные связи в матрицах 
спектральных плотностей, определены элементы главных и побочных диагоналей, учитываются 
корреляционные связи между воздействиями. Эти особенности значительно усложняют вычисления. Поэтому, 
чтобы получить матрицы некоррелированных воздействий, матрица внешних воздействий приведена к 
диагональному виду. 

Результаты исследования. Применением различных методик расчета проведено численное сравнение 
спектральных плотностей и интенсивности колебаний мобильной машины при варьировании скоростей 
движения и характера микропрофиля агрофона. Определенные характеристики спектральных плотностей и 
колебаний мобильных машин агропромышленного комплекса позволили выработать рекомендации по 
практическому применению представленных зависимостей для проектирования данной техники. 

Обсуждение и заключения. Результаты решения матрицы спектральных плотностей внешних воздействий 
различными методами представлены на графике спектральных скоростей колебаний. Анализ графических 
зависимостей показал, что идентичные результаты, вне зависимости от метода расчета, получаются только для 
машин со слабыми функциональными связями при некоррелированном внешнем воздействии агрофона. Для 
некоторых случаев установлены резонансные скорости движения машины. Графически показано влияние 
неровностей агрофона на интенсивность колебаний узлов машины и рассеивание нагрузок на агрегаты. 


Ключевые слова: многомерная динамическая система, матрица спектральной плотности, некоррелированные 
внешние воздействия, обратная матрица, прикрепленная матрица. 
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оп Феш 4пуе ипиз. Уезвик оЁОЗТО, 2020, уо1. 20, по. 2, рр. 155—161. № рз://401.0г2/10.23947/1992-5980-2020-20- 
2-155-161 


Введение. При проектировании надежной мобильной техники необходимо учитывать не только 
статические, но и динамические нагрузки. На характеристики динамических нагрузок значительное влияние 
оказывают колебания машины и ее узлов. Учет и определение колебательных параметров требуют разработки 
специализированных динамических моделей, рассматривающих эти параметры в аспекте вероятностного 
анализа. Основная причина возникновения колебательных воздействий обусловлена неровностью профиля 
пути мобильной машины и непостоянством физико-механических характеристик агрофона. С достаточной 
точностью задачи, учитывающие эти особенности, решаются линейно, что недопустимо при решении 
многомерных динамических систем с большим количеством взаимных связей. 

Материалы и методы. На качество технологического процесса, выполняемого мобильной 
сельскохозяйственной машиной, влияют как характеристики ее рабочих органов [1-5], так и характеристики 
ходовой системы [6,7] и микрорельеф агрофона [8, 9]. В общем случае внешние воздействия на машину 
рассматриваются как система с и функциями. Процессы изменения микронеровностей пути (или изменения 
плотности почвы) обозначим как функции 41 (Е) ...9 (Е) [9, 10]. 


Если статистическое описание внешних воздействий допускает стохастическое ортогональное 
представление, то справедливо уравнение: 


ч() =т [4.0] + [0 (9)-$. (вв) 45, 


где О;(@) — система случайных функцийс; ф; (1,0) — система неслучайных функций { и параметра 02; т [4;(1)] 
— математическое ожидание. 
Мобильные машины АПК — это системы с запаздывающими воздействиями [2, 6, 9], которые 


характеризуются многими параметрами неровности микрорельефа поля или твердости почвы [2, 8, 9]. Эти 
параметры являются однотипными и могут быть сдвинуты во времени 10. Помимо этого, конструкция машины 
(колеса переднего и заднего моста) влияет на взаимосвязь между воздействиями. С учетом этих особенностей и 
предположением, что математическое ожидание внешних воздействий равно нулю, матрица спектральных 
плотностей внешних воздействий описывается следующим образом [11, 12]: 
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5, (©) 5, (©) р Ку (1®) 
’ (<, 15) = 5, (®)е”“* 5, (©) Сл, (16) 2 (1) 
Саб уя (15) Со (1®) и (<) 


где с. =65с», =е ^%; с =1;с, =е/^* — для переднего ведущего моста; с; =е 


—Лоц . 
5 


= =е”“%. 6, =| —ДИЯ 
заднего ведущего моста. 

Для упрощения вычисления статистических характеристик приведем матрицу к диагональному виду, 
чтобы получить матрицу некоррелированных внешних воздействий. Введем характеристическую матрицу и 
приведем матрицу (1) к виду: 

5, (<) г. ^ (75) 5, (®)е ^^ Сб, (1®) 
+ + Ро, . . 
5, (в, /0)-^(15)-Е=| 5, (®)е^*  5,(®)-^(1ю) сб, (1) |, (2) 
сб (76) — с,5„(1) 5, (6)-№(до) 
где 4 — характеристическая функция; Е — единичная матрица. 
Характеристическое уравнение представим в виде: 


а? (16) +Ь^. (1%) +с^.(1ю)+4=0. (3) 
Характеристические функции ^. ( 1®),^., (15), ^. ( 1®) определялись с использованием формул Кардано 

[13]. 
Матрица перехода к новому базису, построенная по координатам векторов ^. ( 1®) (при 2 =1, 2, 3), 


определенных из уравнения (3), имеет вид: 








Ти 1 ТП 
Т(/2)= "> п» Пы|- (4) 
Тз 13 Пз 
В этом случае обратная матрица имеет вид: 
ч( и) 70°) 
р (15) —- “(6 > (5) 
где Г(/%) — присоединенная матрица; |Т ( 9) — определитель матрицы. 
Искомая матрица определялась по выражению: 
5, (18) =Т ' (1®)-5, (в, /2)-Т (15) = 41а |5, (16)-5> (16) 5, (15) , (6) 








где 51, 52, 93 — спектральные плотности некоррелированных воздействий функций комплексного переменного. 
Нахождение диагональной матрицы (6) упрощается решением плоской задачи. В этом случае внешние 
воздействия представим в виде: 














| 5, (©) 5, (®)е 
5, (в, /ю)=| и, } 7 
9 (о ©) К (=) ее“ в. (©) ( ) 
Тогда характеристическая матрица определится выражением: 
5, (®)-^(°) 5, (®)е"* 
5 —^(<).Е=|" ь | 8 
Ч (с, 5) (<) о (©) е®% 8. (©) —^.(®) ( ) 














В свою очередь, характеристическое уравнение представим как: 

^^ (©)-25, (®)-^=0. (9) 
Откуда корни уравнения равны: /, (©) = 0,^,, (©) =25, (©) : 
При корне 4, (®)=0, система уравнений для определения координат собственных векторов имеет 


следующий вид: 
5, (®) пи, (12) +5, (®)е ^^ "Пр? (1°)=0. (10) 


ЛЮ 


Задавшись 1 ‚ находим 1|> (16) =1. 


— е 
При корне 4. (©) = 25, (©) ‚ получим систему уравнений: 
| (©) 1, (762) +5, (в) е`^* п» (25) = 
5 (в) пы (1) +5, (®)-т (76) =0 


Тогда матрица перехода к новому базису принимает следующий Вид: 


(О 
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— 64 Го 
е 71° е/ о 


1 (22) = (12) 


—1 е?/®% 


Преобразованием выражения (8) найдем спектральные плотности некоррелированных воздействий: 
. | — 16 — 
5 (1) = 5, (®)-(+е ^^) =0. 
ил) — Л) — 
5, (15) = 5, (©). 1+е^*)=0. (13) 
Если допустить, что профиль агрофона в продольном и поперечном сечении определен 
некоррелированными случайными процессами [8,9], то можно представить третий элемент диагональной 


матрицы, как 5, (©) . В этом случае матрица (13) примет вид: 
5, (16) = ав |5, (© )-(1+е ^^); 5, (©). @+е^*); 5, (<}] (14) 


Решением матрицы (1), совместно с матрицей (14) и с учетом формул, описывающих некорреляционные 
воздействия, определялись спектральные плотности колебаний. В итоге получилось выражение: 


5. (©) = 5, (в) [Ф!, (©)+Ф› ()|=0. (15) 

Результаты исследования. На рис. | представлены графики спектральных плотностей вертикальных 

скоростей передней части молотилки комбайна семейства «Вектор». Очевидно, что динамические параметры 

машины влияют на спектр колебаний (рис. 1). Узкополосный спектр характерен для условий близкого 

совпадения частот собственных колебаний мостов комбайна с жаткой (рис. 1). Также на ширину спектра влияет 
запаздывание воздействий от колес мостов комбайна. 


$(6), -103 (м2/с3) 
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1 1,2 1,4 1,6 1,8 2 2,2 2,4 2,6 2,8 3 
©, -10 (с!) 
Рис. 1. Спектральные скорости колебаний при решении с использованием: 
1 — матрица (1); 2 — матрица (6); 3 — матрица (15) 


Результаты расчета интенсивности колебаний с использованием матриц (1), (6), (15) графически 
представлены на рис. 2, что подтверждает зависимость функциональных связей в машине от характера внешних 
воздействий. 
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Рис. 2. Графики интенсивности колебаний при движении машины 
по микропрофилю с различным характером неровностей: 
1 и 2 — вдоль борозды; 3 — поперек борозды поля 


Партко С. А. и др. Проектирование мобильных машин моделированием динамических нагрузок на узлах их приводов 


Рано 5. А. её а. Море тасйте 4е п @йтоизй аупатис [оа4 ятшайоп оп шей! апуе ипИ$ 





Решение систем, представляющих собой некоррелированные внешние воздействия на движущийся 
машины со слабыми функциональными связями, дают близкие результаты вне зависимости от используемых 
методов. 

Одинаковые результаты при использовании различных методов расчетов получаются только для машин 
со слабыми функциональными связями, при движении по микропрофилю поля, создающего некоррелированное 
внешнее воздействие. 

Оценка динамических характеристик комбайнов производилось по следующей формуле: 


И, = 15, (6-4, (16) 
По 


где 1, 2. — границы частот спектра. 
Интенсивность колебаний для различных условий, определенная с использованием формулы (16), 
представлена на рис. 3. 
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Рис. 3. Результаты инженерной оценки динамики комбайна 


На рис. 3 представлены зависимости перегрузок — от скорости движения комбайна. Очевидно, что 


различный спектральный состав неровностей микропрофиля агрофона влияет на интенсивность колебаний 
узлов машины. Следует отметить наличие резонансных скоростей движения в интервале от 10 до 15 км/ч. 
Полученные графические зависимости (рис. 2, 3) позволяют оценить разброс нагрузок на агрегаты. Затененная 
область иллюстрирует их рассеивание и позволяет оценить коэффициенты вариации (2). В таблице приводятся 
параметры распределения нагрузок при средних скоростях движения машин [2]. 
Таблица 1 
Динамические параметры нагрузок для мобильных машин АПК 























Тип Скорость 2. 

Режим работы ы 
машины движения м/с 8 

й 4 
Комбайн «Вектор» И - — __ 
транспортный 2,781 0,835 
МТА трактор К-700А, тележка рабочий 3,352 0,474 
2 ПТС-4 транспортный 8,323 0,510 
МТА трактор рабочий 1,247 0,364 
Т-150А, комбайн ККУ транспортный 3,527 0,382 
МТА Т-150А 

и. 29 Рабочий (пашня) 0,752 0,346 




















Обсуждение результатов. Анализ результатов расчетов спектральных плотностей колебаний по 
формулам (1), (6), (15) графически представлен на рис. 2, 3. Из графических зависимостей видно, что спектр 
воздействий нагрузок и их динамика оказывают влияние на погрешность вычислений. 
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Расчет динамических многомерных систем и оценка колебаний машин автоматизирован и производится 
с использованием специализированного программного продукта’. Программа применяется для выбора 
колебательных параметров контуров и агрегатов мобильных сельскохозяйственных машин; расчета плавности 
хода зерноуборочных комбайнов; для проектирования деталей узлов машин [14]; оценки условий труда 
механизатора. Перед использованием программного продукта необходимо предварительно выбрать 
динамическую модель; задаться конструктивными параметрами машины и ее агрегатов; определить массово- 
геометрические и упруго-диссипативные характеристики системы. 

Выводы. Выражения (1), (6), (15) приемлемы для проектировочных расчетов спектральных плотностей 
воздействий на узлы машины. Использование формулы (6) недопустимо при проектировании машин с режимом 
движения по неровностям агрофона с низкочастотным спектральным составом и для машин с низкой 
корреляцией колебаний ходовой и управляемой части. 
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Влияние массы деталей на удельный съем при виброабразивной 
обработке 
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Введение. Следует указать на недостаточность исследования проблемы влияния массы деталей на процесс 
виброабразивной обработки. В трудах А. П. Бабичева и М. А. Тамаркина отмечается факт такого воздействия, 
но не раскрываются степень и механизм влияния. В формулах съема металла учитывается только число 
взаимодействий, приводящих к микрорезанию. Цель представленной работы — определить влияние массы 
деталей на удельный съем металла при виброабразивной обработке. 

Материалы и методы. Использован эмпирический, т.е. экспериментальный подход. В качестве образцов 
выбраны детали из материалов Д16 и ЗОХГСА, широко применяющихся в авиационной промышленности. Для 
изменения массы в заготовках просверлены отверстия, в некоторые образцы залит свинец, в другие забиты 
заглушки, выполненные из того же материала, что и сами заготовки. Таким образом, эксперименты 
проводились с образцами цельными, полыми, а также утяжеленными свинцом. Рабочая абразивная среда — бой 
шлифовальных кругов 40х80 мм, зернистостью 25 и трехгранные призмы 15х15 мм, зернистостью 16. 
Эксперименты позволили наглядно показать влияние зернистости среды на удельный съем обрабатываемой 
детали. 

Результаты исследования. Определены параметры влияния массы деталей на удельный съем при 
виброабразивной обработке. Полученные результаты позволяют показать съем на единицу площади. Данные 
аппроксимированы методом наименьших квадратов линейной функцией. Вариант ее распределения выбран с 
помощью статистического критерия Фишера. 

Обсуждение и заключение. Показано, каким образом при виброабразивной обработке масса обрабатываемых 
деталей обусловливает удельный съем. В дальнейшем следует пополнять базу данных, которые используются 
при определении влияния характеристик обрабатываемого материала на рассматриваемый процесс. Это 
позволит ввести поправочный коэффициент влияния массы в формулу съема металла, что обеспечит более 
точное прогнозирование съема металла на стадии проектирования технологических процессов виброабразивной 
обработки. 


Ключевые слова: виброабразивная обработка, абразивная среда, отделочно-зачистная обработка, масса, 
удельный съем, аппроксимация, критерий Фишера, коэффициент влияния массы, бой абразивных кругов, 
зернистость. 
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Введение. Перед машиностроением всегда стояла задача повышения качества выпускаемой продукции. 
С развитием технологий совершенствуются способы ее решения. Для этого необходимы методики 
прогнозирования эффективности процессов финишной обработки. Виброабразивная обработка позволяет 
обеспечить требуемые параметры качества при высокой производительности, обрабатывать детали сложной 
формы, а также множество деталей одновременно. С целью решения рассматриваемых вопросов изучаются: 

— организация и развитие технологических процессов и методов воздействия на рабочую абразивную 
среду и предмет, подлежащий обработке или переработке; 

— создание новых сред и способов обработки материалов; 

— снижение энергозатрат; 

— улучшение качества обработки [1]. 

При виброабразивной обработке металл и его окислы снимаются с поверхности вследствие взаимных 
соударений частиц среды и заготовок. Этот процесс обеспечивается вибрацией рабочей камеры, в которой 
находятся заготовки и среда. Камера установлена на пружинных опорах, поэтому может колебаться в 
различных направлениях. Колебания сообщаются от инерционного (или другого вида) вибратора с частотой до 
50-100 Гц и амплитудой от 0,5-5,0 мм и более [2]. Количество взаимодействий, приходящихся на единицу 
поверхности детали в единицу времени, носит случайный характер [3-10]. 

Цель исследования — определить влияние массы деталей на удельный съем металла при 
виброабразивной обработке. 

Материалы и методы. Обработка проводилась на универсальном вибрационно-галтовочном станке с 
четырьмя рабочими камерами объемом 10 литров. Для замеров массы использовались аналитические весы АД 200. 

Рабочей абразивной средой был бой шлифовальных кругов 40х80, зернистостью 25 [11] (рис. 1), а также 
призмы трехгранные (ПТ) 15х15 мм, зернистостью 16 (рис. 2). 





Рис. 1. Бой кругов 40х80, зернистость 25 
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Рис. 2. Призмы трехгранные, зернистость 16 


После черновой токарной операции образцы были обработаны в среде трехгранных призм в течение 
10 минут, чтобы удалить заусенцы и выровнять шероховатость (рис. 3). 





Рис. 3. Обрабатываемые детали из алюминия Д16 и стали ЗОГСА 


Затем производилась обработка в два этапа по тридцать минут в среде боя абразивных кругов. Рабочая 
камера вибрировала с частотой 34,7 Гц и колебалась с амплитудой 2,5 мм при непрерывной подаче 
технологической жидкости (раствор кальцинированной соды 0,2 %). Раствор удалял продукты износа (частицы 
металла и абразива) с поверхности деталей и рабочей среды. Затем проводилась обработка в среде ПТ (также в 
два этапа по тридцать минут). 

Результаты исследования. В результате экспериментов были получены искомые значения массы т, г 
и удельного съема. Оценены отклонения 4 от теоретической модели. Рассмотрены в сравнении расчетные и 
табличные значения критерия Фишера (табл. 1—4). 

Таблица 1 
Результирующие данные показателей массы и удельного съема образцов из стали ЗОХГСА в среде ПТ 






































Масса т, г Удельный съем, г/мм” а (отклонение) 

28,53175 2,18976Е-06 —1,28317Е-07 
28,5348 1,73433Е-06 —5,83836Е-07 
40,60125 2,8835Е-06 2,10629Е-07 
40,6065 2,99742Е-06 3,24395Е-07 
4472755 2,98196Е-06 1,87783Е-07 
4473305 3,12387Е-06 3,2954Е-07 
78,52885 3,29591Е-06 —4,91895Е-07 
78,5414 4,01564Е-06 2,27469Е-07 
96,7192 3,47528Е-06 —8,4725Е-07 
96,48725 4,55972Е-06 2,44012Е-07 
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Масса т, г Удельный съем, г/мм” а (отклонение) 
103,5209 4,14608Е-06 —3,76394Е-07 
103,5378 5,42683Е-06 9,03865Е-07 
Невязка 2,68051Е-12 
Стандартное отклонение 1,01462Е-06 
Доверительный интервал 95 % |9,09679Е-07 
Среднее значение 3,40253Е-06 
Правая граница интервала 4,3122Е-06 
Левая граница интервала 2,49285Е-06 

















Е расчетное 32,24585006 
Етабличное 3,105806516 
Угловой коэффициент а 2,93961Е-08 
Свободный член В 1,47936Е-06 




































































Масса т, г Удельный съем, г/мм а (отклонение) 
28,55685 3,41274Е-06 2,74176Е-08 
28,5631 3,55469Е-06 1,69061Е-07 
40,6304 3,99566Е-06 1,37735Е-08 
40,6367 3,59643Е-06 —3,85769Е-07 
44,7538 4,0894 1Е-06 —9,6225Е-08 
44,7618 4,54399Е-06 3,57958Е-07 
78,58825 6,05373Е-06 1,96296Е-07 
78,6054 5,48683Е-06 —3,71445Е-07 
96,5534 6,59611Е-06 —1,49006Е-07 
96,775 6,5642Е-06 —1,91866Е-07 
103,5847 7,8442Е-06 7,51664Е-07 
103,6058 6,77172Е-06 —3,21859Е-07 
Невязка 1,21985Е-12 
Стандартное отклонение 1,5276Е-06 
Доверительный интервал 95 % |1,3696Е-06 
Среднее значение 5,20914Е-06 
Правая граница интервала 6,57874Е-06 
Левая граница интервала 3,83954Е-06 
Е расчетное 8,4308727 
Е табличное 3,105806516 
Угловой коэффициент а 4,94112Е-08 
Свободный член ВБ 1,97429Е-06 









































Масса т, г Удельный съем, т/мм” а (отклонение) 
10,13095 1,36343Е-06 3,59001Е-07 
10,13305 1,19318Е-06 1,88631Е-07 
14,59555 1,10785Е-06 —1,5278Е-07 
14,59835 1,5057Е-06 2,44909Е-07 
27,98755 1,00601Е-06 —1,02312Е-06 
27,9912 1,1657Е-06 —8,63641Е-07 
34,66113 1,83039Е-06 —5,81706Е-07 
34,73845 4,14834Е-06 1,73181Е-06 
60,9715 3,36435Е-06 —5,57569Е-07 














Таблица 2 


Результирующие данные показателей массы и удельного съема образцов из стали ЗОХГСА в среде боя кругов 


Таблица 3 


Результирующие данные показателей массы и удельного съема образцов из алюминия Д16 в среде ПТ 
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Масса т, г Удельный съем, г/мм” а (отклонение) 

60,98585 4,57721Е-06 6,54468Е-07 
Невязка 6,11719Е-12 

Стандартное отклонение 1,36499Е-06 

Доверительный интервал 95 % |1,40279Е-06 

Среднее значение 2,12622Е-06 

Правая граница интервала 3,529Е-06 

Левая граница интервала 7,2343 1Е-07 

Е расчетное 13,93008274 

Е табличное 3,249835542 

Угловой коэффициент а 5,713851Е-08 

Свободный член Б 4,23063Е-07 








Таблица 4 


Результирующие данные показателей массы и удельного съема образцов из алюминия Д16 в среде боя кругов 


















































Масса т, г Удельный съем, г/мм” а (отклонение) 

10,12275 3,18169Е-06 8,85432Е-07 
10,1236 2,31113Е-06 1,47567Е-08 
14,58795 3,63568Е-06 7,25513Е-07 
14,5936 3,20988Е-06 2,98936Е-07 
27,97185 3,54486Е-06 —1,20542Е-06 
27,9844 4,00792Е-06 —7,44089Е-07 
34,43015 4,09858Е-06 —1,53964Е-06 
34,43015 6,08322Е-06 4,45003Е-07 
60,92935 1,10349Е-05 1,75338Е-06 
60,99307 8,65641Е-06 —6,33871Е-07 
Невязка 9,45129Е-12 

Стандартное отклонение 2,79828Е-06 

Доверительный интервал 95 % |2,87576Е-06 

Среднее значение 4,97643Е-06 

Правая граница интервала 7,85219Е-06 

Левая граница интервала 2,10066Е-06 

Е расчетное 51,65197229 

Е табличное 3,249835542 

Угловой коэффициент а 1,37487Е-07 

Свободный член В 9,04509Е-07 

















Аппроксимируем табличные данные линейной зависимостью методом наименьших квадратов. 
Возьмем приближающую функцию в виде у = ах + Б. Тогда невязка (сумма квадратов отклонений) имеет вид: 
5$(а,Б) = УР (у; - ах, — Ь)?. В методе наименьших квадратов невязка должна быть минимальной. В точке 
минимума функции нескольких переменных частные производные этой функции по независимым параметрам 


равны нулю, следовательно, условия минимума: 


== о —ах— Бух; = 0, 


о, ах -Ь) =0. 


После преобразований получаем следующую систему двух алгебраических уравнений с двумя 
неизвестными: 
ау ей. ЕВ 11 =У- 124: 1 
о (1) 
аут 1% +5т = у 1У. 


Бутенко В. И. и др. Влияние массы деталей на удельный съем при виброабразивной обработке 
Вшепко И 1. е а. Ерес! ор та55 ор рат оп гетоуа пие ипаег уБтоабтаяте тасёття 





Обозначим х значения массы деталей, у — удельный съем деталей. Аппроксимируем заданную 
таблично функцию линейной зависимостью. Для определения наилучших параметров а и Б методом 
наименьших квадратов решим систему (1). Решим систему матричным методом в среде Мисгозой Ехсе| и 
получим значения а и В (см. табл. 1-4). 

Для проверки адекватности результатов воспользуемся критерием Фишера и занесем их в табл. 1- 4. 
Расчетное значение критерия Фишера имеет вид: 


Е * ХО: расчет. — Усреднее ращег п-ЕЁ-1 
асчет. —^ 
р Е ХО - У раечет 





где { — количество факторов х, влияющих на у, и — количество наблюдений. 

Сравнивая расчетное и табличное значения коэффициента Фишера (см. табл. 1-4), видим, что Е 
расчетное существенно превосходит табличное. Таким образом, можно сделать вывод, что с достоверностью 
95 % построенная зависимость соответствует исходным данным. 

Подставим полученные значения а и 6. 

Для образцов ЗОХГСА в среде ПТ у = 2,93961Е — 08 хх + 1,47936Е -— 06, в среде боя кругов у = 
4,94112Е — 08 хх + 1,97429Е -— 06. 

Для образцов Д16 в среде ПТ у = 5,73851Е — 08 хх + 4,23063Е - 07, в среде боя кругов у = 
1,37487Е — 07 хх + 9,04509Е - 07. 


В графическом виде результаты представлены на рис. 4-7. 
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Рис. 4. График зависимости от массы удельного съема материала образцов ЗОХГСА, обработанных в ПТ 
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Рис. 5. График зависимости от массы удельного съема материала образцов ЗОХГСА, обработанных в бое кругов 
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Рис. 6. График зависимости от массы удельного съема материала образцов Д16, обработанных в ПТ 
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Рис. 7. График зависимости от массы удельного съема материала образцов ЗОХГСА, обработанных в бое кругов 


Из графиков видно, что удельный съем значительно меняется в зависимости от среды, а также от массы 
и материала обрабатываемых деталей. 

Заключение. Анализируя рис. 4—7, можно сделать вывод, что при обработке деталей большей массы 
удельный съем увеличивается. Это связано с изменениями количества движения (импульса). Так как импульс 
равен произведению массы тела на его скорость, то при увеличении массы детали увеличивается импульс 
взаимодействия частиц с поверхностью деталей. При увеличении массы деталей в два раза удельный съем 
увеличивается в 1,5-2,2 раза. При сравнении значений углового коэффициента а и свободного члена В между 
обработками можно сделать предположение, что на них влияют значения параметров зернистости сред, а также 
характеристики обрабатываемого материала. Полученные результаты позволяют уточнить теоретические 
зависимости, введя коэффициент для определения влияния массовых соотношений. Это обеспечит более точное 
прогнозирование съема на стадии проектирования технологических процессов виброабразивной обработки. 
Полученные результаты пополняют базу данных, которые используются при определении влияния 
характеристик обрабатываемого материала и среды на рассматриваемый процесс. 
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штамповочного инструмента для различных условий нагружения 


С. Д. Колотиенко, А. В. Журавлев, Е. В. Рощина 


ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 


Введение. Статья посвящена аналитическому определению общих критериев для прогнозирования 
износостойкости тяжелонагруженных пар трения, в частности, штамповочного инструмента, работающих при 
различных условиях нагружения. В основу предложенной системы критериев положены физические 
зависимости, связывающие основной показатель износостойкости — число циклов нагружения Мк с 
параметрами трещиностойкости п и С. 

Целями работы являлись разработка и аналитическое обоснование расчетной методики, позволяющей 
прогнозировать износостойкость штамповочного инструмента для различных условий нагружения, а также 
создание предпосылок для экспериментального подтверждения эффективности применения разработанной 
методики. 

Материалы и методы. Предложены математические модели, связывающие основной критерий 
износостойкости — число циклов нагружения № с параметрами трещиностойкости и и С. Проведено 
аналитическое обоснование предложенных моделей. 

Результаты исследования. Разработаны математические модели для прогнозирования износостойкости 
штамповочного инструмента, работающего в различных условиях нагружения. В частности, для случаев: при 
трении скольжения и качения инструмента о пластически деформируемый металл, для условий 
термомеханической контактной усталости, многоцикловой хрупкой контактной ’повреждаемости, 
пульсирующего контакта о пластически деформируемый металл, а также для комбинированных 
разновидностей фрикционного контакта. 

Обсуждение и заключения. Полученные результаты могут быть использованы в процессе проектирования и 
оптимизации конструкций штамповочного инструмента, работающего в различных условиях нагружения, а 
также при прогнозировании долговечности его работы. 


Ключевые слова: штамповочный инструмент, критерии износостойкости, условия нагружения, контактная 
повреждаемость, трещиностойкость, скольжение, качение, наплавочные материалы. 
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Введение. Экспериментальные и теоретические исследования процессов изнашивания наплавочных 
материалов свидетельствуют, что никакие свойства материала в отдельности не определяют однозначно 
износостойкость в условиях экстремального характера фрикционного контакта. Только определенный комплекс 
физико-механических свойств с учетом особенностей фрикционного взаимодействия должен учитываться при 
определении износостойкости наплавочных материалов. 

Критерии износостойкости должны отвечать условиям конкретных трибосопряжений. Для определения 
критериев следует придерживаться в каждом отдельном случае единого подхода. Очевидно, что при всех 
явлениях и разновидностях контактной повреждаемости следует учитывать главное условие: тот материал 
считается более износостойким, у которого более длительна стадия роста усталостных трещин. Таким образом, 
критерий износостойкости должен учитывать продолжительность процесса разрушения поверхностного слоя и 
отделения первой частицы износа. В этом случае определение критерия износостойкости сводится к 
определению числа циклов №», при котором происходит разрушение поверхностного слоя. 

В этом случае назначение критериев износостойкости может состоять не в расчете контактной 
выносливости материалов, а в сравнении и подборе наплавочных материалов для конкретных условий 
фрикционного контакта, характерных для тяжелонагруженных пар трения. Таким образом, критерии должны 
обеспечивать адекватность рядов износостойкости для конкретных разновидностей фрикционных контактов. 

В соответствии с единым подходом при определении критериев обходимо учитывать следующие 
основные положения. Во всех видах фрикционного контакта следует рассматривать один тип трещины. В 
общем случае обязательно должны учитываться параметры трещиностойкости материалов, которые 
определяют кинетику процессов развития усталостных микротрещин. При сравнении наплавочных материалов 
между собой по критериям обязательно необходимо учитывать возможность развития конкурирующего вида 
износа — малоцикловой вязкой повреждаемости в результате лепесткового отслаивания. 

Материалы и методы. Для условий реверсивного скольжения о пластически деформируемый металл 
автором впервые был предложен метод анализа износостойкости с учетом основных свойств материала, 
которые определяют процесс разрушения поверхностного слоя при трении [1]. 

Оказалось возможным определить критерий износостойкости с учетом вязкости разрушения, предела 
текучести, коэффициента трения, размеров исходных трещин. Как показывают исследования [2-9, 11, 13, 15, 
16], данный подход является перспективным и для других условий фрикционного контакта, характерных для 
тяжелонагруженных пар трения, в том числе штамповочного инструмента. 

В. В. Рубанов, исходя из критерия, предложенного для реверсивного скольжения в работе [1], и 
используя зависимость предела текучести от ряда характеристик структуры материала, предложил определять 
износостойкость материала по критерию [4]: 


5.441166 |-К, 
Г 


К, г а 2 (1) 


где о, предельное механическое напряжение в материале от действия матрицы, при котором происходит 
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его разрушение; /; — средний размер дефекта; К‹— вязкость разрушения; ф — объемная доля дефектов; © — 
модуль сдвига;  — вектор Бюргерса; /— коэффициент трения (здесь и далее применяется коэффициент 
трения на основе закона Зибеля). 


Учитывая, что для ряда сплавов существует зависимость К 1с ОТ Предела прочности при изгибе ос после 


изг 2 
ряда упрощений В. В. Рубанов рекомендует для оценки износостойкости наплавочных материалов критерий: 
С`ф-5„, 
те — Го .(1-$).2? . (2) 
р 

Данный подход имеет ограниченные возможности, так как предложенный критерий является общим 
для различных схем фрикционного контакта и при определении износостойкости наплавочных материалов для 
скольжения и качения по критерию (2) не учитываются различия в механизме контактной повреждаемости, в 
условиях контактного взаимодействия для скольжения и качения с тангенциальным усилием или без него. 

Критерий износостойкости (2) не учитывает параметры трещиностойкости материалов, которые 
непосредственным образом влияют на скорость развития усталостных микротрещин, а, следовательно, 
определяют длительность процессов контактной повреждаемости. 

Результаты исследования. При аналитическом определении критериев износостойкости следует 
принять основное допущение: поле контактных растягивающих напряжений однородно, а его величина равна 
максимальному значению, которое определяется на поверхности. Таким образом, данное допущение не 
учитывает постепенного затухания растягивающих напряжений в глубине материала и является более жесткой 
оценкой. Как показали исследования [10], возникающие при этом погрешности невелики, если длина трещин 
меньше глубины расположения зоны сжимающих напряжений. 

Задача определения критериев износостойкости материалов, изнашивающихся в условиях контактной 
повреждаемости, состоит в следующем. Необходимо аналитически определить длительность (число циклов 
контактного взаимодействия) процесса отделения частицы износа недопустимых размеров. В связи с этим в 
качестве растягивающих напряжений следует использовать те, которые в конкретном фрикционном контакте 
воздействуют на трещину. 

Учитывая различия между малоцикловой хрупкой контактной повреждаемостью и многоцикловой 
повреждаемостью, задача определения критериев износостойкости разбивается на два типа. В первом случае 
для малоцикловой хрупкой контактной повреждаемости при определении длительности процесса контактного 
разрушения следует учитывать возможность достижения трещины критического значения. Во втором случае 
для малоцикловой хрупкой контактной повреждаемости, как установлено, трещины не достигают критических 
значений (/.>>1[,). В этом варианте возможен случай использования в качестве конечной длины трещины 
предельно допустимой величины размера выкрашивания (устанавливается либо из требований допуска на 
износ трибосопряжения, либо из требований ограничения размера контактных повреждений на кромке рабочей 
поверхности). 

Критерий износостойкости для скольжения о пластически деформируемый металл 

Длительность процесса контактного разрушения определяется интегрированием уравнения Пэриса 
[17]. В качестве величины растягивающих напряжений для этого случая фрикционного контакта следует взять 
[10]: 
-4 2 (3) 
ртах д 2 
гдер — среднее давление в контакте. 

В этом случае выражение для коэффициента интенсивности напряжений примет вид: 


ма АР. (4) 


Длительность процесса разрушения, а, следовательно, и критерий износостойкости Мк, определится как 
оС" / 








№» = = ь 5 
Е ] СК""? С(арГ)! 7 ( ) 
При п-2 2 решение выражения (5) примет вид 
2" 11 (6) 
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Из основных положений механики разрушения известно, что критическая длина трещины, это 
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величина, при достижении которой трещина начинает закритическое самопроизвольное движение. Данный 
момент однозначно определяет вязкость разрушения материала К\с. В случае контактного нагружения 
однородным полем напряжений при скольжении значение критической длины трещины может быть 
определено из выражения [10]: 


И" у? 





Е (1) 
После соответствующей замены и преобразований, получим 
2”? т (4 \"" 
(8) 





п-2 
2 


мкл 


1С 


М» = 
(-2)С (ру || 
7 
Данное выражение можно использовать в качестве критерия износостойкости при скольжении о 
пластически деформируемый металл в условиях высоких контактных давлений, когда реализуется 
малоцикловая хрупкая контактная повреждаемость. В этом случае, как следует из выражения (7), 
износостойкость наплавочных материалов определяется контактным давлением, коэффициентом трения, 
вязкостью разрушения, размерами дефектов структуры и параметрами трещиностойкости п и С. 
Критерий износостойкости для термомеханической контактной усталости 
Воздействие на поверхность трения циклов нагрев-охлаждение, характерное прежде всего для 
пульсирующего контакта о нагретый пластически деформируемый металл, приводит к возникновению в 
поверхностном слое периодических термомеханических напряжений. Известно выражение В. И. Духовченко 
[12] для определения этого вида напряжений: 
а счЕАТ | (9) 
1-у 
где о — температурный коэффициент линейного расширения; ЕЁ — модуль упругости; АТ — температурный 
перепад цикла «нагрев-охлаждение»; у — коэффициент Пуассона. 
С учетом выражения (9), решение интеграла для уравнения Пэриса при п: 2 в этом случае будет 
определяться по формуле: 


_ 2А-ъ)" |: 
(п-2)Ст"?(аЕАТ)" = р. 


1 








(10) 


Е 





В случае принятия за [с предельного значения В тах» которое назначается в качестве допуска на 
предельный размер очага выкрашивания кромки рабочей поверхности, получим 
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2 п2 
(п-2)Ст"? (аЕЛТ)"1,? Ь,2. 
Данное выражение может быть использовано в качестве критерия износостойкости наплавочных 
материалов при термомеханической усталости. Износостойкость материалов в этих условиях также зависит от 
параметров трещиностойкости п и С, определяется величиной допуска на единичное выкрашивание, модулем 
упругости, коэффициентом Пуассона и обратно пропорциональна коэффициенту теплового расширения, 
величине температурного перепада и размеру исходных дефектов. 
Критерии износостойкости для многоцикловой хрупкой контактной повреждаемости 
Учитывая в качестве [= [; и принимая во внимание, что при качении с тангенциальным усилием и 
упругом пульсирующем контакте угол развития трещин не зависит от свойств материала и близок к 45°, для 


определения длительности процесса разрушения при качении с тангенциальным усилием получено выражение: 


п 


2 7,35 1 1 
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Для качения без тангенциального усилия 


2 [462 т 








г п- п- 2 (14) 
ев ы. |" 
При пульсирующем контакте 
2 [6,65 | 11 
В п-2 п-2 (15) 
(п = 2)С О ах 12. 12. 


Так как в условиях многоцикловой хрупкой контактной повреждаемости трещины при развитии не 
достигают критического значения, /.>>[;, то выражения (13)—(15) можно преобразовать следующим образом: 
— для качения с тангенциальным усилием 


2.7,35" 
М, = В (16) 


(п-2)С(+ИГ)" о" [2 


р тах 7 





— для качения без тангенциального усилия 











м, = 2.4, 62 — (17) 
("-2) Со" „и? 
— ДЛЯ пульсирующего контакта 
= 2-6, 65 = (18) 
(п —-2) Со ьЁ 2 


Полученные критерии износостойкости свидетельствуют о том, что и в этом случае определяющее 
влияние оказывают параметры трещиностойкости материалов и и С, размер дефектов, а, кроме того, в случае 
качения с тангенциальным усилием, следует учитывать коэффициент трения. 

Критерий износостойкости для качения о пластически деформируемый металл 

Особенность этой разновидности фрикционного контакта состоит в том, что контртело приобретает 
интенсивную пластическую деформацию. В этом случае напряжения в зоне контакта не могут оцениваться в 
соответствии с формулами для упругого контакта и критерии износостойкости (16)—(18) неприменимы. 

Для решения данной задачи воспользуемся аналогией для качения о пластически деформируемый 
металл и скольжения в этих же условиях. Тогда на трещину будет периодически воздействовать пластически 
деформируемое катящееся тело. Учитывая, что качение в этих условиях сопровождается значительным 
проскальзыванием, воспользуемся выражением (6) для контакта скольжения. В качестве конечной длины 
трещины следует использовать величину предельного допуска на единичное выкрашивание кромки рабочей 
поверхности. 

В этом случае число циклов для разрушения поверхностного слоя можно определить через выражение: 








Е (19) 
(п-2)С(4рт) [2 Ь2 
После преобразований выражение (19) можно переписать следующим образом: 
2” ьт - Г? 
м, = (20) 


(п-2)С(4р/)" 1? В 
Таким образом, износостойкость наплавочных материалов при качении о пластически деформируемый 
металл определяется параметрами трещиностойкости и и С, контактным давлением, коэффициентом трения, 
предельным допуском на размер очага выкрашивания и размерами исходных дефектов. 
Критерий износостойкости для пульсирующего контакта о пластически деформируемый металл 
Характерными чертами этого вида фрикционного контакта являются чрезвычайно высокие контактные 
давления (до 4000 МПа), что связано с динамическим характером приложения контактной нагрузки. В этом 
случае также невозможно воспользоваться выражением (18) для упругого пульсирующего контакта. 
За максимальные растягивающие напряжения, действующие на трещину примем величину, которую 
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при первом приближении можно определить исходя из условия о, „„=р по формуле: 


1-2 
бк = г р (21) 





гдер — средняя величина контактного давления. 
Длительность процесса разрушения определяется интегрированием уравнения Пэриса для данных 
условий контактного нагружения: 
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Представим значение критической длины трещины Г через Ж. : 
К’ 
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3 
и после преобразований получим: 
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1 

Данное выражение представляет собой критерий износостойкости наплавочных материалов для 
пульсирующего контакта о пластически деформируемый металл. В этих условиях износостойкость 
определяется вязкостью разрушения, параметрами трещиностойкости, коэффициентом Пуассона, контактным 
давлением и размерами начальных дефектов. 

Критерии износостойкости для комбинированных разновидностей фрикционного контакта 

В условиях трения скольжения о нагретый пластически деформируемый металл рабочая поверхность 
инструмента подвергается поочередному воздействию напряжений от сил трения и термомеханических 
напряжений. В соответствии с гипотезой о линейном суммировании напряжений коэффициент интенсивности 
напряжений с учетом выражений (3) и (4) в этом случае можно представить в виде: 


4р’ аЕАТ 
К=ои м = В т (25) 
и. 1-у 
Тогда интегрированием уравнения Пэриса можно получить выражение для определения необходимого 
для разрушения числа циклов: 


у 11 2 1 1 
М; =| д аЕАТ п — д аЕЛТ п п-2 п-2 | (26) 

) НЕ: 

"С р п п? (п 2)С М, 12 [2 

п 1-у п 1-у 
Критическую длину трещины можно определить из выражения (25) имея в виду, что момент 
закритического разрушения наступит при К = в. . Тогда 
К/ 
Ра ео. (27) 
«(427 „оват) 
И 1-у 











Подставляя значение /‹ в уравнение (26), можно получить выражение для критерия износостойкости в 
случае комбинированного воздействия на поверхность трения скольжения и термомеханических напряжений: 











п-2 п-2 
| п 2 (Е +алАт 
ИЯ —у 
-- 4рГ аЕАТ\"| “= КГе-Ъ 28) 
п"? (п | р Е [2 1с 
—у 


Комбинированное воздействие на поверхность трения возможно при упругом пульсирующем контакте 
с наложением термомеханических напряжений. В этом случае уравнение (18) следует записать в виде: 


п 








2 6, 65 1 1 
= | (29) 
Е п-2 п-2 
(п-2)С а и аа 
== 7 Г: 


Учитывая, что в этом варианте нагружения >>», получим выражение для критерия износостойкости 
при упругом пульсирующем контакте с термомеханической усталостью 
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2.6,65" 
п п * 
ви. + маг 


2 

1 
Обсуждение и заключения. Получены аналитические зависимости для расчета показателя 
износостойкости Ме для различных условий работы штамповочного инструмента на основе общих критериев, 
описывающих параметры трещиностойкости материалов п и С: Полученные результаты могут быть 


М, = 





(30) 








1-у 


использованы в процессе проектирования и оптимизации конструкций штамповочного инструмента, 
работающего в различных условиях нагружения, а также при прогнозировании долговечности его работы. 
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Введение. Для повышения оперативности принятия решений на предприятии целесообразно использовать спе- 
циальный программный комплекс интеллектуальной поддержки. Такой продукт необходим при проектирова- 
нии системы защиты информации и повышении ее неуязвимости в ходе модернизации или изменения конфигу- 
рации. Цели исследования: создание алгоритма и математической модели программного комплекса интеллек- 
туальной поддержки принятия решений. 

Материалы и методы. Метод поддержки принятия решений при проектировании системы защиты информации 
базируется на использовании нейронной сети (многослойный персептрон). Для объективной оценки исходной 
защищенности информационной системы (ИС) сформирована математическая модель анализа событий без- 
опасности. 

Результаты исследования. Проанализирована статистика злоумышленных воздействий на ИС предприятий. 
Определена необходимость своевременной и точной модернизации системы защиты информации. Важными 
характеристиками процесса проектирования системы защиты информации являются скорость получения ре- 
зультата и снижение остаточного риска ИС. В связи с этим актуально использование систем искусственного 
интеллекта в процессе определения лучшего набора подсистем защиты. Классифицированы угрозы нарушения 
информационной безопасности (ИБ). Определены основные классы событий безопасности. Создана математи- 
ческая модель нейронной сети, указаны входные параметры ее работы. 

Действующая ИС предприятия генерирует многочисленные события, что обусловливает необходимость авто- 
матического сбора и анализа данных с подсистем регистрации объектов ИС. Детально рассмотрен процесс ана- 
лиза событий безопасности, так как от корректности данных, полученных таким образом, зависит адекватность 
сгенерированных проектных решений. Сформирован алгоритм работы программного комплекса. 

Обсуждение и заключение. Полученные результаты могут быть использованы при проектировании системы 
защиты информации на предприятии. Кроме того, администраторы ИБ могут применить разрабатываемый 
программный комплекс для корректировки конфигурационных настроек средств защиты информации. 
Предложенное решение позволит минимизировать деструктивное влияние разработчика системы защиты, 
который может быть и бывает субъективен. 


Ключевые слова: информационная безопасность, информационная система, нейронная сеть проектирования, 
многослойный персептрон, алгоритм. 
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а|2отИт ап4 а табетайса! по4е] оЁ Ве зоЙ\’аге расКасе ог ниеПеста[ 4ес1$10п зиррог. 

Маепаб апа Мейо4$з. ТВе 4ес11юп заррой тефо4 ипаег 4е51е01т$ ап шоппаНоп зесигИу зузет 1$ Базе оп Ше 
и5е оРа пеига| пебмогК (тиШауег регсерёгоп). Рог ап об]есйуе аззеззтеп{ оЁ Ве ша зесигИу оЁ ап шЮгтаНоп 
зумет (15), а па етайса| плоде! ог Ше апа1у$1$ ОР зесигИу еуеп 15 4еуе]ореч. 

Везий5. ТВе зайзис$ оЁ таНс1ои$ аНасК$ оп Ше 1$ оРГещегривез 15 апа[уте4. ТВе пее4 Юг Ите]у ап4 ассигае тод- 
егпт7айоп оЁ Ве шЮгтайоп ргоесйоп зузет 15 дейегитеа. Ппроцапе сВагасепзИс$ оЁ Пе 4ез1епте ап шЮппайоп 
зесигИу зузет аге е зрее4 аё \\сП {Пе гези 15 оМаше4 апа е тедисноп ш е гез1ма[ 115К оЁ 1$. ш 1$ гегага, 
фе изе оГагийсла! пцеШзепсе зузетз ш е ргосез$ оЁ деегите фе Без зе оЁргойесйоп зибзузетлз 15 ппрогаге. 
ТБе {йгеаё5 ю субег зесигиу (СЗ) аге Сазуе4. ТВе таш с1аз$ез оЁ зесигИу сует аге 4ейпе4. А тафетайса! то4е1 
оЁ Фе пеига] пе\уогК 15 4еуе]оре4; ше шри{ рагатеег$ ог Из орегайоп аге ш@саеа. Тве ситгепе егиегризе 1$ сепег- 
аёез питегои$ еуепё$ \БлсЬ песеззЦа(ез фе ашютайс соЦесНоп ап4 апа[у$1$ оР даа Нот забзузетлз Юг гео15егте 
15 оБ]есё5. ТБе ргосез$ оЁ апайу7ште зесигИу еуеп5 15 сопз14еге ш де зшсе Ше адедиасу оЁ Пе сепегае4 4е51еп 
Чес11015 4ереп@$ оп бе сотестез$ оЁВе даёа оате4 ш 15 \’ау. ТВе а!еогИбт оЁ Пе зоЁ\аге расКазе 1$ Когтед. 
Обсизяюоп апа Сопсиз1опз. ТВе гези сап Бе изе4 ш е 4ез1еп оЁ Фе шЕЮгтаНоп зесигИу зузет аё Фе етмегризе. 
ш аа4юп, С$ адтиизгаюг$ сап зе Фе 4еуе!оре зоЁ\уаге расКазе №0 а4л15ё Ше сопЯгигайопт зе 12$ оЁ 
шЮгтаноп зесигйу 10015. ТЬе ргорозе4 зойоп \Ш плшитхе Фе 4езгасйуе шИчепсе оЁ Фе 4еуеорег оЁ те 


зесигИу зузет \/шсН тпау ап4 Варреп 1ю Бе заБ]есйуе. 
Кеупога$: субег зесигиу, поппайоп зузет, пеига| пебмотК епошееппе, тиауег регсербоп, а1еогибт. 


Еог сйапоп: Е. А. Упепбтге, А. У. МПазроуа. А1еогИт оЁ зоЙ\уаге расКаее оЁ ниеПесвла! 4ес1$1юп зиррог \Пеп 
Чезюпшя субег зесигиу зузбет аё Фе ещегризе. Уезник оЁ ОЗТО, 2020, уо|. 20, по. 2, рр. 178—187. 
Б&рз://401.0го/10.23947/1992-5980-2020-20-2-178-187 


Еипатя тргтапоп: Фе гезеатсЬ 1$ 4опе \ИБ Фе Йпапс1а| заррог оЁ уоцпе Киззап $с1еп55 Нот бе Огап 
Соипсй оф Ргез1Аеп оЁ 1е Киззап Еедеганоп ш В &0)О оп “ВиЙате а тоде! оЁ пиеПесла| 4ес11оп-такте зирроге 


уУВеп 4ез12птФ ап шРоппайоп зесигИу зузет аё бе ещегризе”. 


Введение. С развитием промышленности в России растет число предприятий, отнесенных к объектам 
критической информационной инфраструктуры (ОКИИ). Статистика распространения атак от вендора РозШуе 
Теспо!оз1ез' показывает рост количества успешно реализованных злоумышленных воздействий в этой сфере. 
В 2019 году зафиксировано 125 атак на промышленные информационные системы (ИС). Это более чем в три 
раза (или на 212%) превосходит аналогичный показатель 2018 года (40 атак). Диаграмма распределения коли- 


чества атак по кварталам года (О) приведена на рис. 1. 


'Актуальные киберугрозы: ТУ квартал 2019 года / Розуе Тесвпооз1ез. ОВГ.: ВИрз://Ам\иу.р!зесигИу.сот/та-га/гезеагсЬ/апа1унс/субегзесшу- 
{геазсаре-2019-44 (дата обращения: 24.02.2020). 
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Рис. 1. Число атак на ИС предприятий в 2018 и 2019 году 


Промышленные ИС атакуют преимущественно с использованием вредоносного программного обеспе- 
чения (90 % атак). Данную сферу курируют Федеральная служба безопасности и Федеральная служба по техни- 
ческому и экспертному контролю. Необходимость обеспечения информационной безопасности (ИБ) промыш- 
ленных информационных систем подтверждается статистикой ведущих аналитических центров [1]. Согласно 
нормативно-правовым актам Российской Федерации’, владельцы ОКИИ задействуют комплекс организацион- 
но-технических мероприятий для обеспечения безопасного функционирования информационной инфраструк- 
туры. При этом законодательством предусмотрена периодическая ревизия качества функционирования системы 
защиты информации (СЗИ), оценка ее эффективности. Этим обусловлена необходимость оперативной коррек- 
тировки конфигурационных настроек имеющихся средств или доукомплектация системы инструментами защи- 
ты информации. Особую роль играют оперативность и точность принимаемых решений, а величина остаточно- 
го риска не должна превышать установленные показатели [2]. В связи с этим предлагается автоматизировать 
процедуру поддержки принятия решений при проектировании СЗИ на предприятии. 

Материалы и методы. Для поддержки принятия решений при проектировании СЗИ применен метод, 
основанный на нейронной сети (многослойный персептрон) [3]. Входные данные для работы нейронный сети 
— угрозы нарушения информационной безопасности и события безопасности. Кроме того, для объективной 
оценки исходной защищенности информационной системы сформирована математическая модель анализа со- 
бытий безопасности. В рамках этой модели рассчитанные меры сходства ИС сравниваются с одним из уровней 
безопасности ИС. В качестве метрики сходства используется взвешенное расстояние Манхэттена. 

Результаты исследования. В рамках данного исследования предложен подход к проектированию си- 
стемы защиты, основанный на показателях важности подсистем защиты информации, входящих в СЗИ [4]. С 
учетом этих данных предлагается формировать перечень средств защиты наиболее важных подсистем. Данный 
подход позволяет усилить СЗИ за счет нейтрализации актуальных угроз. При этом актуальность следует опре- 
делять, исходя из автоматизированного анализа событий безопасности [5, 6]. Вектор важности подсистем защи- 
ты информации формируется в соответствии с выражением 

У = 5(С1а55_Тйт), (1) 
где И = (И,,..., 5) — вектор важности подсистем СЗИ; С1а5$_Тйт = (С1а55_Ту, ..., С а5$_ТИть) — вектор ак- 
туальности классов угроз; 5 — функциональная зависимость, определенная нейронной сетью. 


2 О безопасности критической информационной инфраструктуры Российской Федерации : федер. закон от 26 июля 2017 г. № 187-ФЗ / Гос- 
ударственная Дума РФ, Совет Федерации // ФСТЭК России. ОВГ: Б@рз://Ещес.галекисвезКауа-хазбсВИа-шогпай/обезресвеше- 
Бегоразпози-киисвезКо]-шЁогла1юоппо]-птйазиаКгу/285-хаКопу/1610-Ее4егашту]-гаКоп-о{-26-1ушуа-2017-2-п-187-Ё2 (дата обращения: 
18.05.20). 
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Множество угроз ТАт нарушения ИБ для удобства разделим на классы С1а55_ТИт, которые определя- 
ются, как в [6]: 

С1а5$_Тйт = { Трт_Вт, ТАт_ Г, ТА_О15Е, ТАг_Го55$, Трт_В, ТАт_А}, (2) 
где Трг_Вт — класс угроз нарушения ИБ типа «взлом»; ТИг_Ё — класс угроз нарушения ИБ типа «утечка»; 
ТП_015Е — класс угроз нарушения ИБ типа «искажение»; ТАт_[.05$ — класс угроз нарушения ИБ типа «утра- 
та»; ТАг_В — класс угроз нарушения ИБ типа «блокирование»; Тйт_А — класс угроз нарушения ИБ типа «зло- 
употребление». 

Для формирования вектора актуальности угроз нарушения ИБ ИС предлагается формировать две мат- 
рицы соответствия: 
— множества угроз множеству классов угроз — МТйг, 


— множества событий безопасности классам угроз — МЕтепё. 


Матрица МТйт: 
МТит = Тйт х С1а55_Тйт = (тё)). (3) 
Здесь тёй;; определяется по формуле: 
пей, = | если угроза принадлежит классу угроз, (4) 
0 в противном случае. 
Здесь п — количество классов угроз, которым принадлежит угроза. 
При этом УЕ, У, тёН;; = 1. 
Матрица МЕуепЕ: 
МЕтетё = Ететё$ Хх СТйт = (тет,,). (5) 
Здесь тет; определяется по формуле: 
1, если событие возникает при реализации угрозы класса, 
Е | 0 в противном случае. (6) 


Для формирования вектора важности подсистем защиты информации предлагается использовать в про- 


граммном комплексе нейронную сеть — многослойный персептрон, функционирующий согласно формуле, 
приведенной в [7]: 
ок — Йк 
Оиё; = ГОлил Пил — 6%), (7) 


Пул = бий п 
где ок — К-Й нейрон входного слоя; 7, — К-й элемент входного вектора; Оиё;; — выходное значение /-го 
нейрона 1-го слоя; { — функция активации нейрона, которая определяется функциональной зависимостью 
(1); ил — вес [го входа /-го нейрона Г-го слоя; 1; — значение [-го входа /-го нейрона 1-го слоя; 9;, — уро- 
вень активации /-го нейрона 1-го слоя; Оиё;_1, — выходное значение /-го нейрона (1-— 1)-го слоя. 

На этапе анализа событий безопасности исходными данными являются журналы событий безопасно- 
сти, которые создает системное и прикладное программное обеспечение ИС предприятия. Множество событий 
безопасности Еуеи{5 включает в себя следующие классы [8]: 

Еуетё5 = {ЕщетЕт, Мападетептё5 иБЕт , Ассез5ОБ}Еу ,РоПсуСвпапдеЕу, 
ОзеРтищедезЕт, [5Ртосе$5езЕт, Гете 5Ез}. (8) 
Здесь ЕтетЕт — события класса «вход субъектов в систему»; Мападетете5 иБЕу — события класса «управ- 
ление субъектами»; Ассез5ОБ]Ет? — события класса «получение доступа к объектам»; РойсуСпапдеЕт — со- 
бытия класса «изменения политики системы»; И5еРти?ЦедезЕт — события класса «использование субъектом 
особых привилегий»); [5Ртосе55езЁт — события класса «функционирование процессов системы»; Гете $Ет 
— события класса «уровень системы». 

Из множества событий выбираются опасные. Для сопоставления множества опасных событий с клас- 
сами угроз используются матрицы соответствия. Матрица, сформированная в рамках данного исследования, 
представлена в табл. 1. 
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Таблица 1 
Матрица соответствия угроз нарушения ИБ ИС предприятия 
Типы угроз 
ее 
а 5 ъ 
3 |5 3 ыы 
я 1 о Е г 5 
ь Е 5 ы ы - — 
Класс события безопасности ра 5 а 5 © о а 
Ф = |” ‚2 > в > 
з о 1% о 5 ъ ъ 
> 5 % > 2 © 5 
ы З в а [м + к 
> А = 
5 < сз в 9% 
з о «а — 
= 5. > 
ТАт_Вт + + + 
Тйи_ Г, + + 
ТИ_215 + 
ТА"_Г05$ + + 
Тит_В 
Тик _ А + + + + 
































Затем формируется вектор актуальных классов угроз нарушения ИБ ИС. После этого актуальным клас- 
сам угроз противопоставляются подсистемы СЗИ. Соответствие определяется с помощью нейронной сети с 
учетом функциональной зависимости У в соответствии с формулой (1). 

Немаловажным фактором является количество анализируемых наборов событий и их источников. При 
этом стоит отметить, что количество событий безопасности прямо пропорционально количеству источников — 
информационных ресурсов ИС предприятия [8]. С учетом большого числа событий, генерируемых работающей 
ИС, актуально проводить автоматический сбор и анализ данных с подсистем регистрации объектов ИС, описы- 
вающих события. В связи с этим стоит подробно рассмотреть анализ событий безопасности, так как от кор- 
ректности полученных на этом этапе выводов зависит адекватность сгенерированных проектных решений. 

Согласно [9, 10], в определенный момент времени Т текущее состояние ИС предприятия 5т Е 
5тате, 5тайе = {5потт, 5аапд, 5апотт} можно охарактеризовать как: 

— нормальное (5п07т)) — штатное функционирование системы в соответствии с ее задачами и соглас- 
но документам, регламентирующим работу; 

— опасное (54ап.д) — некорректное функционирование ИС, фиксируются нарушения работы, связан- 
ные с атаками злоумышленника, сбоями и отказами программного и (или) технического обеспечения; 

— аномальное (бапотт) — временное изменение штатного режима функционирования ИС и всплеск 
аномальной активности пользователей, программ и сетевого трафика. 

Более детального анализа требуют опасные и аномальные события. Набор такого рода событий — это 
входные данные этапа сопоставления классов угроз. Такие события свидетельствуют о реализации угроз без- 
опасности. 

Следовательно, входными данными для нейронной сети (7) будут результаты мониторинга и анализа 
событий безопасности. 

Любое событие безопасности ЁуепИ5, может быть описано кортежем атрибутов [9]: 

Еуептё[$; = (ТО, Бава, Гете|, боитсе, ЕуетЕТуре, Етепё5 Ее, беситеРатат5), (9) 
где 2 — код события; ра — время генерации события; Геуе! — уровень опасности события; Зоигсе — ис- 
точник возникновения события; Еует!Туре — тип события; Еуеп5 кие — состояние события; ЗесигеРатату — 
вектор параметров безопасности события. 


5еситеРататз = (1, и, 715К), (10) 
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где й — показатель возникновения события определенного кода относительно общего количества событий без- 
опасности за период ДТ; и — тяжесть последствий события (потенциальный ущерб); 715 — риск нарушения 
информационной безопасности ИС. 

За период АТ в информационной системе генерируется множество событий ИС (ЕуепИ5), которые 
необходимо оценивать для определения уровня защищенности ИС. Соотношение числа событий того или иного 


типа к общему количеству определяет показатель р: 


в = 


МЕтеп О 
“МЕуети Но 
где МЕуеп/ Ш — количество событий определенного кода, МЁуеи! — общее число событий за период ДТ. 
Прикладные программные продукты по-разному определяю коды событий. Так, в рамках данной рабо- 
ты рассмотрена кодировка ОС \/тдо\из корпорации Мисгозой?. 
Риск нарушения информационной безопасности ИС — это функция от частоты реализации события и 


потенциального ущерба: 


715К = Пхи. (12) 
Сумма частных показателей позволяет определить общий риск А 15Кбит: 
ВК = УТ К, (13) 


Множество, обеспечивающее классификацию событий и оценку частных показателей безопасности 
ИС, описывается подмножествами элементов: 

РР5 = {{ЕуТуре}, {Еубшие}, {155ие}}. (14) 
Здесь ЕуГуре — множество типов обнаруженных событий, Еу5кие — множество возможных состояний собы- 
тий: 

Еубкие = {Еут, ЕуЧ, Буа}, (15) 
где Еу" — события нормального функционирования ИС; Е?“ — события нарушения ИБ ИС; Ет“— аномаль- 
ные события, характеризующие отклонения ИС от штатного режима функционирования (нуждаются в допол- 
нительном анализе). 

Для определения принадлежности события ЁуепИ5; к одному из трех состояний Еу5 ие решается зада- 
ча классификации — используется множество 5Р = РатдеЕу У МоттЕъ, разделяющееся на два базовых под- 
множества. Они формируются экспертной группой на основании данных о типовых событиях штатного режима 
функционирования ИС и ранее обнаруженных атаках и инцидентах ИБ, описывающих сигнатуру типовых для 
данной ИС и опасных событий: 

— РатдеЁу — множество событий, которые являются известными признаками атаки или определяют сцена- 
рий развития инцидента; 
— МоттЕу — множество событий, характерных для штатного режима функционирования ИС. 
Е(Етепт 5$; (Ет5 саде) > {5Р};) классифицируются следующим образом: 
Ет", Етет Г 5; Е МоттЕъ, 


Етет 15; (Ет5таде) = ЕуЧ, ЕтетЕ 1$; Е РатадеЕт, (16) 
Ез“, Еует 15; & МоттЕу 9 РатадеЕт. 


Здесь ЕтетЕ 15; (Еу5таёе) — события ИС с атрибутом состояния, каждому из которых соответствует набор свя- 
зей 5Р; с типовыми событиями из множества шаблонов штатного режима и режима с нарушением ИБ. 

Если событие отсутствует в профилях штатных событий или событий нарушения безопасности, оно 
определяется как аномальное. Причины его возникновения должны отдельно рассматриваться администрато- 
ром системы ИБ. 

155 ие ={15погт, 5аапе,15апогт} — множество состояний ИС, указывающее на штатный режим функ- 


ционирования ИС. Состояние ИС определяется по формуле: 


ЗОписание событий системы безопасности в У шао\$ 7 и \Ушдо\з Зегуег 2008 В2 / Мисгозой. ОВГ: Вирз://зиррог.писгозой.сот/го- 
гае!р/977519/4езсирноп-оЁЕ-зесигИу-еуеп-ш-\уи140%5-7-ап4д-ш-улио\з-зегуег-2008 (дата обращения: 05.08.2019). 
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5потт, УЕтет 15; Е Еуепт!5 |Ет5 ще = Етп, 


155а4е = 3 [54апд, ЗЕуеп $; Е Еуеп 15 |Еубъаке; = Бу, (17) 
5апотт, [5потт п 15аапд, 


где Гбиогт — режим нормального функционирования ИС; /5апогт — режим функционирования ИС с призна- 


ками аномальной активности; /54апе — режим с зафиксированными нарушениями безопасности ИС (рис. 2). 


Описывается Описывается 
множество множество 
нормальных 1$дапо событий 
событий ТЗапопи безопасности 
Описывается 
множество 
аномальных 
событий 


Рис. 2. Отношения подмножеств состояний ИС 


С учетом полученных данных формируется вектор, определяющий состояние ИС, суммарный риск, а 
также долю аномальных событий и событий нарушения безопасности: 

15 = (155 ае, Е15Кзит, МапЕт, МаапаЕъ), (18) 
где Мапа[Ет — доля обнаруженных аномальных событий, М4апаЁт — доля событий нарушения ИБ ИС пред- 
приятия. 

Для принятия решения о необходимости усиления подсистем защиты информации рассчитывается уро- 
вень безопасности ИС 5Г: 

51, = {безопасное, стабильное, аномальное, кризисное, опасное}. 

Для определения уровня безопасности ИС проектировщик системы ИБ формирует вектор, определяю- 
щий эталонный показатель защищенности ИС 15_Рет}. 5Ё формируется на основании сходства векторов 15 и 
15_Рет}. 

Для оценки безопасности ИС в рамках данной математической модели принято пять уровней, поэтому 
при определении принадлежности ИС к одному из уровней применен метод А ближайших соседей [10]. С этой 
целью составляются вспомогательные векторы, соответствующие оставшимся четырем уровням безопасности. 
Выполнение данной операции базируется на значениях вектора 1[5_Рет}{, соответствующего безопасному уров- 
ню состояния ИС при умножении на скалярный корректирующий коэффициент (значение определяется экс- 
пертным путем): 

Кх В1$ = (Кх 155аще, К х В15Кзит, К х МапщЕтете, к хх МаапаЕтети). (19) 

В качестве метрики сходства используется взвешенное расстояние Манхэттена [10]: 

р(15, 15) = и У.4_ 115; — В!$). (20) 

На основании представленной математической модели разработан обобщенный алгоритм работы про- 


граммного прототипа (рис. 3). 
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Рис. 3. Алгоритм работы программного прототипа комплекса интеллектуальной поддержки принятия решений при 
проектировании системы защиты информации на предприятии 


Шаг 1 (блоки 1, 2). Запуск программного прототипа. Начало алгоритма. Ввод входных данных о векто- 
рах требований 15_Ретг} к состоянию ИС каждого уровня безопасности. Ввод периода, через который будет 
проводиться мониторинг. Загрузка из БД шаблонов с множествами нормальных и опасных событий. 

Шаг 2 (блоки 3—5). Мониторинг журналов событий ОС, получение записей о каждом событии. Форми- 
рование списка событий. 
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Шаг 3 (блок 6). Анализ данных о собранных событиях, формирование кортежа событий — формула (8). 
Классификация событий по формулам (14—15) на нормальные, аномальные и опасные. Расчет частоты возник- 
новения событий, потенциального ущерба и риска по формулам (11-13). 
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Шаг 4 (блоки 7-11). Классификация состояний ИС по формуле (17) на основе данных о распределении 
событий на множество аномальных, нормальных и опасных (шаг 3). 

Шаг 5 (блоки 12—14). На основании рассчитанных данных и кортежей событий в ИС формирование 
вектора текущего состояния, вычисление мер сходства между векторами /[5 и [5_Рет{: формулы (19), (20). При- 
нятие решения о принадлежности ИС к одному из пяти уровней безопасности. 

Шаг 6 (блок 15). Завершение алгоритма. 

Заключение. В рамках данного исследования предложен подход к моделированию системы защиты 
информации. Учитывается различный качественный и количественный состав средств защиты в зависимости от 
актуальных угроз нарушения информационной безопасности. Представленный метод позволяет повысить эф- 
фективность работы внедряемой системы защиты информации и уменьшить вероятность ошибки проектиров- 
щика. Программный комплекс, создаваемый в рамках данного исследования, обладает преимуществами, кото- 
рые невозможно получить при «ручном» проектировании: 

— учет всех данных о защищаемой системе, 
— получение точных результатов в максимально короткий срок. 

Библиографический список 

1. Майорова, Е. В. Методические аспекты реагирования на инциденты информационной безопасности 
в условиях цифровой экономики / Е. В. Майорова // Петербургский экономический журнал : [сайт]. — 2020. — 
№1. — ОВГ: Б&рз://субейепшКа.го/атие/ллеюсвезке-азреКбу-геазтоуатуа-па-ш1Чещу-шЕЮюгта1юппоу- 
Беторазпоз$Н-у-изоутуаН-#5тоуоу-еКопош! К (дата обращения: 24.02.2020). 

2. Применение методов теории нечетких множеств к оценке рисков нарушения критически важных 
свойств защищаемых ресурсов автоматизированных систем управления/ А. И. Братченко, И. В. Бутусов, 
А. М. Кобелян, А. А. Романов // Вопросы кибербезопасности : [сайт]. — 2019. — №1(29). — ОВС: 
Б&рз://субейеттКа га/атие/л/ргилепеше-теюдоу-{еоги-песве-тпо7Незу-К-о{зепКе-пзКоу-пагазвешуа- 
КиисвезК1-уа7впуЮ-зуоуз6у-тазс1зсваетуН-гезигзоу (дата обращения: 24.04.2020). 

3. Витенбург, Е. А. Математическая модель интеллектуальной поддержки принятия решений при про- 
ектировании системы защиты информации на предприятии / Е. А. Витенбург // Промышленные АСУ и кон- 
троллеры. — Москва : Научтехлитиздат, 2019. — С. 54—60. 

4. Витенбург, Е. А. Выбор элементов комплекса защиты информационной системы предприятия на 
основе требований нормативно-правовых документов / Е. А. Витенбург, А. А. Левцова // Вестник Донского 
государственного технического университета. — 2018. — №3. — С. 333-338. 

5. Степанова, Е. С. Разработка модели угроз на основе построения нечеткой когнитивной карты для 
численной оценки риска нарушения информационной безопасности / Е. С. Степанова, И. В. Машкина, 
В. И. Васильев // Известия ЮФУ. Технические науки: [сайт]. — 2010. — №11. — С. 31-40. — ЧВГ: 
Б#рз://субе|епшка га/агас1е/п/гатгабоКа-то4еП-иэго7-па-озпоуе-розгоетуа-песпеоу-Коступоу-Кагу-Фуа- 
сы$еппоу-о{зепК1-п5Ка-пагазвептуа-шогтпа1оппоу (дата обращения: 24.04.2020). 

6. Витенбург, Е. А. Модель угроз информационной системы предприятия/ Е. А. Витенбург, 
А. А. Левцова // Промышленные АСУ и контроллеры. — 2018. — №9. — С. 46-50. 

7. Бова, В. В. Применение искусственных нейронных сетей для коллективного решения интеллекту- 
альных задач / В. В. Бова, А. Н. Дуккарт// Известия ЮФУ. Технические науки. — 2012. — №7 (132). — 
С. 131-138. 

8. Смоляк, Д. С. Мониторинг событий информационной безопасности техногенных объектов / 
Д. С. Смоляк, Т. А. Пулко // Доклады Белорусского государственного университета информатики и радиоэлек- 
троники. — Минск : Изд-во Белорус. гос. ун-та информатики и радиоэлектроники, 2015. — С. 43—47. 








9. Использование методов системного анализа для решения проблемы обеспечения безопасности со- 
временных информационных систем / И. В. Машкина, А. Ю. Сенцова, М. Н. Гузаиров, В. Е. Кладов // Известия 
ЮФУ. Технические науки. 2011. — № 12 (125). — С. 25-35. 

10. Астрахов, А. В. Противодействие компьютерным атакам. Технологические основы. Электронное 
учебное издание : [сайт] / А. В. Астрахов, С. М. Климов, М. П. Сычев. — Москва : МГТУ имени Н. Э. Баумана, 
2013. — 70 с. — ОВГ: Бир://у\у\уса 1 бтазеа.гиЛта 1 О/сотар-айаКк-есбпо.раР (дата обращения: 18.05.2020). 


Сдана в редакцию 27.03.2020 
Запланирована в номер 27.04.2020 


Витенбург Е. А. и др. Алгоритм работы программного комплекса интеллектуальной поддержки принятия решений 


Гйеприга Е. А. её. Аоогийт ор зоДгаге расказе о 1теПесиа! дес1я1оп зиррот зйеп 4еявтия суБег 5есигйу зуйет а йе етегрие 





Об авторах: 

Витенбург Екатерина Александровна, аспирант кафедры «Информационная безопасность», ФГАОУ 
ВО «Волгоградский государственный университет» (400062, РФ, г. Волгоград, пр. Университетский, 100), 
Кезеагспе ШО: М- 0О-8740-2017, ЗсоризШ 57209346586, ОВСШ: |@рз://огс19.0г2/0000-0002-1534-8865, 
е.уцепбиго@ес-гз.ги 











Никишова Арина Валерьевна, доцент кафедры «Информационная безопасность», ФГАОУ ВО «Волго- 
градский государственный университет» (400062, РФ, г. Волгоград, пр. Университетский, 100), кандидат техниче- 
ских наук, Везеагсве ШО: М№-3217-2016, Зсори$ Ш 57201358403, ОВСТО: Вирз://огс14.0г2/0000-0002-0919-2593, 

п 1звоуа.агта@)уо15и.ги 








Заявленный вклад соавторов 

Е. А. Витенбург — формирование основной концепции, определение целей и задач исследования, со- 
здание математического аппарата обработки событий безопасности, формулирование выводов. А. В. Никишова 
— научное руководство, анализ результатов исследования, разработка математического аппарата интеллекту- 
альной поддержки принятия решений, корректировка выводов. 


Все авторы прочитали и одобрили окончательный вариант рукописи. 


Информатика, вычислительная техника и управление 


— 
со 
1 


Бйр://уезииК.Аоп$а.га 


188 


Вестник Донского государственного технического университета. 2020. Т. 20, № 2. С. 188-195. 155М 1992-5980 е155М 1992-6006 
Резий оГ оп 5ие Тесйшс Отуегуйу. 2020. Ио. 20, по. 2, рр. 188—195. 155М 1992-5980 е155М 1992-6006 





ИНФОРМАТИКА, ВЫЧИСЛИТЕЛЬНАЯ 
ТЕХНИКА И УПРАВЛЕНИЕ 
МГОВМАТТОМ ТЕСНМОГОСУ, СОМРОТЕКВ 
ЭСТЕМСЕ, АХО МАМАСЕМЕМТ 


УДК 62-50 В р5://401.0г2/10.23947/1992-5980-2020-20-2-188-195 








Интеллектуальная система мониторинга и управления 
техническим состоянием мехатронных технологических объектов 





А. К. Тугенгольд', Е.А. Лукьянов”, Р.Н. Волошин*, В. Ф. Бонилья* 

+ ?.3 ФГБОУ ВО «Донской государственный технический университет» (г. Ростов-на-Дону, Российская Федерация) 
4 Технологический университет (г. Кито, Республика Эквадор) 

Введение. Исследуются, в основном, цифровые системы, которые управляют обслуживанием отдельных 
технологических мехатронных объектов (ТМО) и комплексов технологических машин. Отмечена 
многочисленность вопросов поддержания надежности состояния и возникающих неисправностей, а также 
многофакторность применения существующих систем мониторинга и диагностики. В связи с этим указано на 
актуальность задач разработки методов и способов обслуживания технологического оборудования, 
позволяющих принимать решения в условиях неопределенности и ограниченности данных. 

Материалы и методы. Для анализа критичности технического состояния сформирована оценка качества 
эффективности автономного управления состоянием устройств. При этом используется метод нейронечеткой 
системы определения совокупного критерия критичности. Предлагается применять этот подход для выработки 
рекомендаций по оснащению технологического объекта необходимыми средствами поддержания общей 
работоспособности и надежности. 

Результаты исследования. Решение позволяет прогнозировать развитие состояния мехатронного 
технологического оборудования, оповещать персонал об аварийных и иных опасных состояниях, при 
необходимости проводить доработку или корректировку управляющих программ. Предусматривается 
выполнение части функций обслуживания технического состояния самим мехатронным объектом, 
т. е. самообслуживание оборудования. Сформулировано понятие «автономное управление техническим 
состоянием». Рассмотрены структура системы и функции управления. Отмечено, что внедрение 
рассматриваемых систем позволит в значительной степени повысить эффективность использования 
оборудования. Результативность автономного управления устройства или ТМО в целом оценивается в 
соответствии со стандартом [$0 13381-1:2004. На основании этого стандарта и данных, представленных ранее, 
построена структура нейронной сети для оценки автономности управления состоянием. Эффективность 
обучения системы рассмотрена с учетом среднеквадратического отклонения выходов сети от целевых значений 
обучающей выборки. 

Обсуждение и заключение. Представлен перечень основных функций управления при разных уровнях 
автономности технического обслуживания: от сигнализации для предупреждения отказа до полной 


автономности технического обслуживания без непосредственного участия оператора. 


Ключевые слова: цифровые системы, автономное обслуживание, управление техническим состоянием, 
критичность состояния. 
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Введение. Уже более 50 лет в области технического обслуживания машин проводятся серьезные 
разработки, которые фокусируются на методах и стандартах поддержания работоспособного состояния 
оборудования. Создаются новые устройства и методы для определения аварийных состояний, идентификации 
неисправностей и прогнозирования остаточного ресурса. Контроль и мониторинг состояния имеют особое 
значение для технологических мехатронных объектов промышленных предприятий — таких, как: станки с 
числовым программным управлением (ЧПУ), машины кузнечного и литейного производства, 
автоматизированные гибкие линии ипр. Методы обслуживания оборудования приобретают различные 
формы [1, 2], у каждой из которых есть плюсы и минусы, влияющие на уровень контроля технологического 
состояния оборудования. Такой контроль и обеспечение надежности машин являются актуальными научными 
задачами. Данной теме посвящены многочисленные исследования. В качестве примера можно привести работы 
[3-14]. Российские и зарубежные ученые рассматривают, в основном, вопросы создания автоматизированных 
систем мониторинга и диагностики. При использовании для автоматизации методов искусственного интеллекта 
создаются процессорные модули, получившие обобщенное название «системы электронного обслуживания (е- 
тайцепапсе)». 

Цифровые системы электронного обслуживания реализуются совокупностью программно-аппаратных 
компонентов, которые предназначены для получения сенсорной информации и оценки состояния оборудования 
в режиме реального времени. С их помощью прогнозируют изменения состояния технологического 
оборудования, его узлов, формируют сообщения о необходимости определенных действий со стороны 
специалистов, в некоторых случаях изменяют параметры работы оборудования. Рассматриваемые системы 
применяются для персональной диагностики и мониторинга технического состояния отдельных мехатронных 
объектов. Кроме того, их задействуют в системе обслуживания группы технологических машин. 
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В фокусе актуальных исследований — создание программ и алгоритмов проведения обслуживания и 
анализа его результатов. На это нацелены разработки в области качества обслуживания и эффективности 
использования оборудования. На современных производствах для сервисного обслуживания применяются 
такие системы, как Саз1р, Ететла, [сазате, ВетоеПаа5епипе|, Пиегтог, ПУТО, 1РО5$, Ргоеиз, МВРО$, \/$ПЕ, 
Тела и др. [15—18]. 

Следует отметить многочисленность вопросов поддержания надежности состояния и возникающих 
неисправностей, а также многофакторность применения существующих систем мониторинга и диагностики. В 
связи с этим остаются актуальными задачи разработки методов и способов обслуживания технологического 
оборудования, позволяющих принимать решения в условиях неопределенности и ограниченности данных. 

В последние годы изменился основной функционал систем контроля, мониторинга и диагностики, 
средств, методов и технологий обслуживания, прогнозирования и анализа состояния оборудования. Они 
переориентированы с анализа технического состояния на предупреждение и обнаружение неисправностей на 
ранних стадиях эксплуатации. Разрабатываются такие способы контроля, которые позволяют минимизировать 
влияние возникающих неисправностей на работоспособность и время простоя оборудования, на количество 
бракованных деталей и частоту отказов. В конечном итоге новые подходы обеспечивают снижение 
экономических потерь. 

Актуальным стал вопрос о разработке нового решения, основанного на интеллектуальности системы. 
Речь идет о создании автоматизированного автономного управления обслуживанием и, как следствие, об 
управлении автономным техническим состоянием оборудования. 

Автономное техническое обслуживание (АТО) [19] — это комплексное определение. Применительно к 
станкам оно подразумевает такое автоматическое выполнение функций сервиса и обслуживания, при котором 
оборудование самостоятельно поддерживает свою работоспособность (самообслуживание). Если же в 
обслуживание необходимо вмешаться человеку, функции сервиса остаются в руках оператора, то есть 
проблемы устраняются оперативно, без привлечения внешних ремонтных служб и до окончания 
производственного цикла. 

Автономность обслуживания технологических мехатронных объектов (ТМО) в ключе, представленном 
ранее, — пока не реализованное, но перспективное решение. Ожидается, что методы искусственного 
интеллекта в сфере диагностики и сервисного обслуживания получат уровень знаний, равный уровню 
обслуживающего персонала. 

Материалы и методы. На современном этапе развития технологий рассмотрим понятие АТО ТМО в 
следующей интерпретации. Цель решения — поддержание или восстановление работоспособности на основе 
оценки состояния компонента и ТМО в целом. Для этого самостоятельное автоматизированное управление 
целенаправленно отображает состояние и необходимые воздействия на компоненты ТМО [20]. 

Как отмечалось ранее, основным подходом для создания и развития автономности при управлении 
работоспособностью технологического оборудования является применение методов искусственного интеллекта 
и интеллектуального управления в условиях неопределенности и нечеткости предоставленных знаний [21]. 

Результаты исследования. В условиях неопределенности рассмотрим процедуру управления как 
зависимость реального состояния от целевого (исправная работа): 

В* =А- В, 
где Аи В — экспертные заключения о реальном и целевом состоянии ТМО. 

Представление результата нечеткого вывода У зависит от: 

— наблюдаемых данных Х, полученных системой диагностики; 

— системы знаний А -> В; 

— наличия коэффициента К, отвечающего за возможную погрешность измерений. 

Следовательно, вывод У принимает вид: 

У=ХхА*+К =ХХх (А >В) +К, 
где У = {У., У,, ..., Уи} — параметры, отвечающие за управление состоянием оборудования и за вывод команд 
для оператора; Х' = {Х\,Х», ..., Х„} — данные о входном состоянии оборудования. 

Используя такую трактовку автономности технического обслуживания, рассмотрим структуру и 
функции АТО ТМО (рис. 1). 
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Рис. 1. Структура системы управления техническим состоянием оборудования 


Здесь сигналы с датчиков от диагностической системы поступают на устройство управления не 
напрямую, а через блок управления техническим состоянием узлов (БУ ТС). 

Узлы блока управления получают данные с центрального узла автономного управления 
технологической системы (ЦУ АУТС), который: 

— координирует работу и обслуживание всех узлов станка; 

— обрабатывает данные, полученные от диагностических систем; 

— на основе этих данных анализирует состояние, прогнозирует развитие неисправностей и принимает 
решения о необходимых действиях по обеспечению бесперебойной работы оборудования. 

В предлагаемом подходе к обобщенным функциям управления техническим состоянием ТМО 
рассматриваются способы управления при помощи сигнализации о возможности наступления отказов, об 
аварийных ситуациях, о необходимости нечеткой полосы работоспособности устройства и др. [22]. 

Электронная система обслуживания (например, е-Мшта МасЬте) может обеспечивать информационную 
связь АУТС объекта через ВС с подразделениями обслуживания и ремонта оборудования и управляющими 
структурами предприятия [23]. 

К функциям автономного управления относится принятие решений по поведению системы при 
неисправности, поэтому система должна обучаться во время работы. Таким образом, методы диагностики и 
мониторинга совершенствуются, что позволит устранять неизвестные ранее неисправности. 

Аналитическими методами сложно (или невозможно) адекватно представить состояния и изменения 
ТМО. Это обусловлено перечисленными ниже особенностями технического обслуживания такого 
оборудования: 

— сложность структуры мехатронного объекта, состоящего из большого числа узлов; 

— многообразие факторов, влияющих на надежность и работоспособность, точность и 
производительность машин; 

— большой диапазон свойств материалов, из которых изготовлены детали; 

— разнообразие требований к функционалу, параметрическим характеристикам и качеству работы 
блоков, модулей и устройств в соответствии с их назначением, точностными характеристиками, реализуемыми 
производственными задачами. 

Концепция е-МтаМасВте обеспечивает управление с учетом данных о реальном состоянии объекта и 
истории эксплуатации. При этом прогнозируется развитие неисправностей для каждого узла. Предлагается 
использование гамма-процентного ресурса надежной работы, если доступны статистические данные или 
применяется пограничная полоса работоспособности [22—24]. 

При разработке АУТС технологических объектов можно предусмотреть ранжирование уровней 
автономности в зависимости от применяемого на предприятии обслуживающего аппарата [25, 26]. Это зависит 
также от критичности машины в общем составе производства. В данном случае критичность определяет, 
насколько важна для производства та или иная его единица. Чем выше критичность, тем больше 
экономических, физических или временных потерь понесет предприятие от простоя этой единицы. Такие 
данные основываются на статистике использования оборудования, числа отказов, времени простоя, стоимости 


запасных частей и изготавливаемых деталей, количества рабочих смен ит. п. 
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При анализе критичности также важно учитывать будущие состояния машины, степень развития 
неисправностей, планирование приобретения запасных частей. Это необходимо для оценки остаточного 
ресурса и расчета времени до возможного отказа в условиях неопределенности и стохастичности данных о 
состоянии оборудования. Оценка критичности позволит выяснить, насколько важным для предприятия 
является данная единица оборудования, какими затратами обернется остановка и восстановление работы. 

Для анализа максимального количества параметров, влияющих на критичность отказов, построим 
систему обобщенной оценки критичности технического состояния устройств ТМО. Для выходного совокупного 
критерия критичности используем оценку на базе нейронечеткой системы АМЕ[ в среде МАТГАВ. Как 
пример рассмотрим критичность отказов шпиндельного узла вертикально-фрезерного станка НААЗ 
ЗММ (рис. 2). 


Ж1 = 0.5 х2 = 0.5 х3 = 0.5 х4 = 0.5 х5 = 0.5 
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Рис. 2. Оценка критичности отказов. Красная вертикальная линия — показатель 0,5: производительность и точность 
оборудования (а); стоимость ремонта оборудования и запасных частей (6); стоимость простоя оборудования (6); количество 
отказов оборудования (узла) за определенный период (г); время на идентификацию и устранение неисправности (0); 

У = | — полное соответствие выходным параметрам (е) 


Как видно из рис. 2, входные воздействия в нейронечеткой системе отвечают за все основные аспекты, 
которых касается работа оборудования на предприятии [27]. В результате на выходе из системы получаем 
обобщенный критерий критичности У, величина которого может варьироваться от0 до 1. Здесь У=0 — 
входные параметры не соответствуют требованиям производства. По данным мониторинга, параметры вышли 
за пределы 0<У<1, но эксперт может сделать заключение о пригодности оборудования (например, при 
У= 0,75). Используя такую методику оценки критичности, при определении работоспособности однотипного 
оборудования можно обозначить диапазон реагирования по возникающим неисправностям и создать новую 
систему принятия решений в условиях неопределенности. Таким образом, наиболее важное (критичное) 
оборудование получает максимальный отклик со стороны оператора и обслуживающего персонала для 
поддержания его работоспособности. 

Методы нечеткой логики позволяют упростить интеллектуальную систему управления, оперативно 
оценивать состояние объекта. Речь идет об автономном управлении техническим состоянием объекта, поэтому 
необходимо разобраться с тем, насколько глубокой может быть автономность. Станок рассматривается вместе с 
оператором, а значит, оборудование управляется автоматически и действиями оператора. Уровни управления 
зависят от функций, возложенных на человека. 

Уровень 1. Контроль. Определяется наличием или отсутствием аварийной неисправности, которая 
приводит к останову оборудования. 

Уровень 2. Диагностика. Проводится анализ данных о состоянии оборудования встроенными или 
навесными сенсорами. 

Уровень 3. Мониторинг. Определяется текущий уровень работоспособности, принимаются 


оперативные или тактические решения. 
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Уровень 4. Оператор выполняет все основные функции без средств автоматизации станка. Они могут 
быть выведены на экран в качестве инструкции. 

Уровень 5. Оператор выполняет часть функций, не требующих отклонения от производственных задач. 

Уровень 6. Автономность сервиса ТМО без непосредственного участия оператора. 

Если необходимо поддерживать работоспособность оборудования на высоком уровне, целесообразно 
задействовать мониторинг (3-й уровень). Использование на 4-м уровне некоторой автоматизации 
представляется переходом к созданию интеллектуальных систем обслуживания. 

Внедрение рассматриваемых систем позволит в значительной степени повысить эффективность 
оборудования. Продуктивность автономного управления устройства или ТМО в целом оценивается в 
соответствии со стандартом [$0 13381-1:2004. На основании этого стандарта и данных, представленных ранее, 


построим структуру нейронной сети оценки автономности управления состоянием (рис. 3). 
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Рис. 3. Структура нейронной сети оценки автономности управления состоянием! 


Эффективность обучения оценивалась с учетом среднеквадратичной ошибки отклонений выходов сети 
от целевых значений обучающей выборки [20]. 

Обсуждение и заключение. В статье представлены основные положения и подходы, позволяющие 
проектировать, создавать и совершенствовать автономные системы управления техническим состоянием 
технологического оборудования. Работа таких систем предполагает значительную неопределенность и 
ограниченность априорных и текущих данных для принятия решений. Этим продиктована необходимость 
внедрения современных технологий и методов искусственного интеллекта. 

Представлен авторский взгляд на роль, функции и построение систем обслуживания ТМО и другого 
оборудования. Предложенный подход основан на использовании цифровых технологий обработки информации 
и принятия решений, что обусловлено ростом сложности ТМО. В каждой технологической системе в той или 
иной мере реализуются принципы открытости, автономности и управления прецизионными движениями. Это 
необходимо учитывать при проектировании процессов обработки деталей и при реализации процессов 
поддержания работоспособности оборудования. 
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Введение. Рассмотрено функционирование искусственного интеллекта (ИИ) в системах хранения данных. 
Определено преимущество его использования при работе с данными как с экономической точки зрения, так и с 
точки зрения безопасности. Целью работы является внедрение искусственного интеллекта в системы хранения 
данных. Основные задачи: описание методов разделения данных, организации их хранения и противодействия 
угрозам безопасности. 

Материалы и методы. Данные, которые необходимо занести на накопители, разбивается на части таким обра- 
зом, чтобы их можно было восстановить, не имея одной из частей. Это необходимо для того, чтобы иметь воз- 
можность доступа и восстановления информации в случае программного или аппаратного сбоя. 

Результаты исследования. Рассмотрена работа искусственного интеллекта при обнаружении угроз безопасно- 
сти. Так как модель подразумевает взаимодействие пользователей с данными, то было выяснено, каким образом 
происходит управление доступом данных, а также приведено описание способа хранения ключей. 

Обсуждение и заключения. Использование искусственного интеллекта при организации хранилища данных 
позволит увеличить скорость работы системы. Искусственный интеллект с встроенными алгоритмами машин- 
ного обучения позволит реагировать на ситуацию, влияющую на состояние системы. Анализ состояния накопи- 
телей позволит избежать возможного аппаратного или программного сбоя. Минимизация человеческого факто- 
ра в функционировании системы способствует улучшению её работы. 


Ключевые слова: искусственный интеллект, пороговое разделение, угроза, машинное обучение, организация 
хранилища, динамическое изменение, ключ, шифрование, резервное копирование, атака, хеш-сумма. 


Образец для цитирования: Жилин, В. В. Искусственный интеллект в системах хранения данных / В. В. Жилин, 
О. А. Сафарьян // Вестник Донского государственного технического университета. — 2020. — Т.20, №2. — 
С. 196-200. 6 рз://401.0г5/10.23947/1992-5980-2020-20-2-196-200 


© Жилин В. В., Сафарьян О. А., 2020 


Агийаа пиеШоепсе т даа $богасе зу$бет5$ 


у.У. 7вййг, О.А. Зааг’уап? 


'Видеппу МИНагу Асадету оЁ Соттитсайоп ($, Рейегзбиго, Виззап Еедеганоп) 
?Роп Зае Тесыиса! Ошуегзйу (Возюу-оп-Ооп, Визчап Еедеганоп) 


тодисйоп. ТВе агайса1 пцеШоепсе (АГ) ре!огпапсе ш 4айа 5югазе зу%етавз 15 сопз1еге4д. \У/Беп мотКте у 
аа, пе адуатазе оЁ 11$ изе Бо ш есопопис {егил$ ап Юг зесигИу 15 д@егттеа. ТВе \уогк оБ]еснуе 15 йе шео4ис- 
оп оР агийсла! иеШзепсе ш айа зюгасе зузетаз. ТВе Кеу (азКз шуо]уе е дезсирНоп оЁ те#о4$ Юг даёа зерага- 
Ноп, огоашхаоп оЁ (5 зюгазе ап соищегасйоп © зесигйу Шгеай5. 

Маепаб апа Мейод$. ТВе даёа фай зво4 Бе Ее4 шю фе 4туе$ 1$ 4гулАе4 шо раг5 зо ай И сап Бе гезюге4 уро 
опе оЁ Ше рагб. ТЫ$ 1$ песеззагу ю Бе а[е №0 ассез$ ап гесоуег шРогтайоп ш йе еуепё оРа зоЁ\аге ог Ваг4\аге 
ГаИге. 

ВезиЙ5. Тве АГ реогтапсе ип4ег д&есйпе зесигИу Шгеай$ 1$ сопз14еге4. Эшсе е то4е| пирНез е п\егасйоп оЁ 
изегз УИ даа, И уаз Копа оп Бо\/ Фе дайа ассез$ сопёго] 1$ сагле4 оц ап Пе Кеуз аге зюге4. 

Р5сизяоп апа Сопсиз1опз5. Те чзе оР АГ т оггап1ите а даа \уагероизе \Ш зрее4 ир е зует. Агайс1а] пиеШ- 
зепсе у БиШ-шт тасбте-еагите а12огИиа$ Ш рго\!4е гезроп те 0 а зйиайоп фа аЁЙес{5 Фе 5е оЁ Ше 5уз- 


Жилин В. В. и др. Искусственный интеллект в системах хранения данных 


ЯрИт И. У. ем. Алая 1теШоепсе т ааа уогазе зубет5 





{ет. Апа[уз1$ оЁ Фе зе оЁ Ше 4пуез \Ш ауо1 а роз е Бага\уаге ог зоЙ\аге гаЙиге. Мшитханоп оЁ фе Ватап 
асюг ш Фе зует орегайоп сопитщез ю Ше пиргоуетеп оЁ $ \огК. 


Кеупога$: агасла1 пиеШзепсе, гезбо!4 зерагайоп, геаф, пасте 1еагипе, зюгазе огоап1хайоп, Чупапис сВапее, 
Кеу, епсгурНоп, БасКчр, абасК, ВазВ. 


Еог сйапоп: У. У. 7, О. А. Заг’уап. Агайсла] пеШоепсе т дайа зюгазе зузеплз. УезииК оГОЗТО, 2020, уо1. 20, 
по. 2, рр. 196-200. В рз://401.ог2/10.23947/1992-5980-2020-20-2-196-200 


Введение. На современном этапе развития информационных и телекоммуникационных технологий че- 
ловек встречается с различного рода информацией. Актуальным является вопрос хранения этих данных в без- 
опасном месте. Хранением информации в зашифрованном виде занимается криптография — наука, изучающая 
способы сокрытия данных и обеспечения их конфиденциальности". 

Во всех реализациях хранения данных информация содержится статично. То есть, данные будут нахо- 
диться там, где они оказались первый раз. На другой диск или другой сектор информация может попасть только 
после удаления еб с первоначального местоположения. Однако данный факт можно выделить как недостаток 
существующих алгоритмов. 

Для повышения эффективности защиты рекомендуется проводить разделение данных [1] с целью хра- 
нения их частей. Это позволяет организовать отказоустойчивые хранилища, в которых организованы алгорит- 
мы восстановления информации в случае какого-либо сбоя, выхода из строя одного из накопителей, а также в 
случае утери одной или нескольких частей данных. 

Если злоумышленник имеет доступ к диску, он сможет получить необходимую информацию. При ста- 
тичном хранении в случае получения доступа к нескольким таким дискам, нелегитимный пользователь может 
воссоздать информацию. Хранение данных в разобщенном виде повышает защищенность информации. Кроме 
того, в больших хранилищах, как правило, используются накопители на жестких магнитных дисках, у которых 
скорость работы (проведения операций чтения и записи) существенно ниже, чем у твердотельных $5) и флэш- 
накопителей °. В связи с этим приходится выбирать между используемым объёмом и скоростью работы. 

Цель работы — описать функционирование искусственного интеллекта в системах хранения данных. В 
работе поставлена задача — описать алгоритмы работы ИИ при выполнении следующих операций: 

® — запись данных на накопители; 

® — обращение пользователей; 

® — анализ хранилища данных; 

» — в случае угроз безопасности информации. 

Влияние искусственного интеллекта на качество работы системы хранения. В современных 
устройствах, таких как смартфоны и компьютеры, разработчики уделяют особое внимание внедрению искус- 
ственного интеллекта. В данных технологиях реализованы алгоритмы машинного обучения, что увеличивает 
скорость работы и сокращает время отклика на проведение часто повторяющейся информации. Следует рас- 
смотреть, каким образом внедрение искусственного интеллекта и машинного обучения в системы хранения 
позволят существенно повысить качество их работы [2]. 

В настоящее время крупные компании и организации начинают внедрять ИИ в свои хранилища дан- 
ных. Широкое внедрение технологий машинного обучения и искусственного интеллекта способствует повыше- 
нию качества работы на уровне управления. Это облегчает работу администраторам сети и хранилища данных 
путем постоянного диагностирования причин перегрузок и снижения трафика, что позволит им заблаговремен- 
но определять потенциально уязвимые сегменты используемой модели. 

Искусственный интеллект подразумевает использование интегрированных алгоритмов глубокого обу- 
чения, которые смогут прогнозировать состояние всей системы и оперативно реагировать на возможные изме- 
нения. Это позволит существенно сократить расходы по ликвидации последствий, вызванных выходом из строя 
оборудования. Кроме того, внедрение искусственного интеллекта в организацию отказоустойчивых хранилищ 
позволит обеспечить их автоматизацию [3]. Под этим подразумевается анализ состояния системы и обработка 
поступающих данных в динамическом режиме. 

Искусственный интеллект и машинное обучение позволят минимизировать вероятность потери дан- 
ных. В совокупности с избыточными массивами независимых дисков такая система увеличивает доступность и 


* Рябко Б. Я., Фионов А. Н. Кринтография в информационном мире. М., 2018. 305 с. ОВГ: Врз://\уилу лесНбоок.п/боок.рНр?14 Боок=1001 
(дата обращения: 04.05.2020). 

? Сравнение $50 и НОР дисков в реальных условиях использования / Хабр. ОВГ: Вирз://абг.соп/ти/розЁ /394135/ (дата обращения: 
04.05.2020). 
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скорость выхода из вынужденного простоя благодаря интеллектуальному восстановлению данных, стратегии 
резервного копирования и переноса необходимых данных [4]. 

Распределение данных по входным параметрам. Для хранения данных в разделенном виде могут 
быть использованы методы порогового разделения. В классических алгоритмах входные параметры — статич- 
ные величины, что является существенным недостатком. Получение доступа к одним данным путём подбора 
входных параметров ставит под угрозу все остальные данные в системе. 

Таким образом, наибольшее предпочтение в рамках информационной безопасности отдается тому ал- 
горитму, который использует различные входные параметры для порогового разделения. Эти параметры могут 
генерироваться случайным образом по определенному алгоритму или зависеть непосредственно от входных 
данных, которые будут соответствующим образом проанализированы для подбора наиболее предпочтительных 
параметров. Генерацией таких параметров будет заниматься непосредственно искусственный интеллект. При 
этом он должен учитывать количество пользователей, которым доступна эта информация. Такие данные как 
параметры порогового разделения и местонахождение первой доли хранятся в базе. Эта база данных находится 
в зашифрованном виде на ключах конкретных пользователей, что соответствует использованию ассиметрично- 
го алгоритма шифрования. 

Параметры предлагаемого алгоритма по обеспечению надежности данных. Искусственный интел- 
лект при распределении данных по накопителям использует следующие параметры: скорость работы накопите- 
лей, их доступность, свободный объем и показатель надёжности. При этом ИИ для вычисления использует па- 
раметры ценности данных, размер и частоту обращений к ним. 

Алгоритм распределения долей данных по накопителям основан на вычислении коэффициентов нако- 
пителей. Для определения скорости работы накопителя 5 необходимо найти среднее арифметическое скорости 
записи 5.ап и скорости чтения 5.т: 

5= 5зап+5чт 
2 

Доступность накопителя А (ауа|ЙабИу) имеет 3 уровня: 0 — недоступен, | — частая ситуация недо- 
ступности накопителя, 2 — редкая ситуация недоступности накопителя, 3 — постоянно доступен. Показатель- 
ный коэффициент накопителя может быть выражен по формуле: 

Кн =5.А.И- В, 
где Г — объём диска; РА (тенаб у) — надёжность, она ранжируется от 1-го до 10-и баллов; А — количество обраще- 
ний к файлу за определённый промежуток времени. 

Коэффициент значимости файла К‹ зависит от уровня ценности (1 — низкая, 2 — средняя, 3 — высо- 
кая). На основании известных атрибутов этот параметр можно вычислить по формуле: 

Кь =5.И:А. 

Рассмотрим способ выбора накопителя хранения остальных файлов. Для определения приоритета 
накопителя используем номер файла тп. в отсортированной таблице по убыванию параметра Къ. Тогда номер 
диска №1, участвующего в выборке для хранения первой доли данных будет найден по формуле: 


а - Пт 
№, = гоипа -;“) то4 п>, 


где 4 — количество дисков, используемых для хранения; Г — количество файлов в системе; п› — количество 
накопителей; гоипа — округление к ближайшему целому. 

Таким же образом определяется диск №, который возможно будет выбран в качестве накопителя для 
хранения первой доли: 

Клтзь = Ку" Катах, 

где Къ — коэффициент разбиваемого файла; Ка.тах — максимальный коэффициент накопителя. 

Далее в качестве №› выбирается тот накопитель, у которого разница коэффициентов Крё.ё И Кн мини- 
мальна. Для хранения первой доли случайным образом выбирается один из найденных дисков №; и №. 

Преимущества рассмотренного алгоритма распределения данных по накопителям. Классически 
диски, объединенные в один массив данных (КА/-массив), работают по статичному алгоритму. То есть при 
обнаружении первой доли данных можно провести процедуру поиска второй доли. Каждый раз, при обнаруже- 
нии очередной доли вероятность определить конкретный алгоритм распределения повышается. В предложен- 
ном алгоритме распределение происходит на основе информационных данных и накопителей. Другими слова- 
ми, задача нахождения всех долей данных является №Р-полной, то есть не решаемой за полиномиальное время. 
Из этого следует, что использование данного алгоритма повышает надежность системы хранения информаци- 
онных данных. 
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Однако, пользователи не создают ключ, он генерируется в автоматическом режиме и хранится на 
устройствах самих пользователей. Таким образом, при попытке расшифровать данные должен использоваться 
конкретный ключ, иначе данная операция завершится неудачно. Местонахождение ключей изменяется динами- 
чески с учетом времени. Такие изменения могут происходить в определенные часы, либо через какой-либо 
промежуток времени [5]. Искусственный интеллект отвечает за местоположение ключа с учетом доступности 
всех устройств, подключенных к системе. Это позволяет усложнить работу злоумышленнику, чьей целью явля- 
ется доступ к информации, хранящейся на накопителях [6]. 

Так как информацией, находящейся в хранилище, оперируют пользователи, следует определять, кому 
из них разрешено проводить операции с данными, другими словами, как разграничивать доступ к данным. 

Сравнение доступа происходит по таблице хеш-сумм '. Если у пользователя есть доступ к данным, то 
они восстанавливаются. В случае отсутствия записи в базе данных о предоставлении прав доступа, субъекту, 
запросившему доступ к объекту, в нём будет отказано. Кроме того, попытка получения доступа будет зафикси- 
рована в лог-файле, а субъект, которому принадлежит хранимая в системе информация, будет проинформиро- 
ван о попытке получения доступа к его данным нелегитимным пользователем или субъектом. Искусственный 
интеллект в данном случае анализирует действия каждого пользователя с целью принятия того или иного ре- 
шения в случае возникновения определенной ситуации. Например, он определяет, является ли запрос ошибоч- 
ным или имеет какую-либо цель [2]. 

На рис. 1 представлен пример предоставления доступа на восстановление данных и отказ в доступе с 
занесением события в лог-файл и уведомлением пользователя. 
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Рис. 1. Предоставление и отказ доступа на восстановление данных 

Рассмотрим функцию управления доступом. Так как в системе основными субъектами являются поль- 
зователи, то вопрос предоставления доступа к данным является особо актуальным. Очень важен полный кон- 
троль над операциями пользователей. Решение о предоставлении доступа к данным принимает распределяющее 
устройство с поддержкой искусственного интеллекта. Использование ИИ в данном случае позволит увеличить 
скорость работы пользователя над его данными. Это достигается путём трансформации наиболее часто исполь- 
зуемой пользователем информации, что позволяет избежать ожидания времени отклика на отправленный за- 
прос о восстановлении данных. 

Рассмотрим динамическое изменение хранилища. В определенных случаях искусственный интеллект 
управляет переносом данных и отвечает за резервное копирование и восстановление данных, в т. ч. с учетом 
времени. При этом данные мигрируют с одного диска на другой, а информация в базе данных также изменяет- 
ся. Это позволяет минимизировать вероятность прогнозирования состояния системы в конкретный момент вре- 
мени. Так как все ключи генерируются на основе предыдущих последовательностей с использованием хеш- 
сумм, задача нахождения ключа, предназначенного для расшифровывания базы данных является МР-полной 
задачей [7]. Искусственный интеллект в данном случае анализирует состояние накопителей, а также системы в 
целом. При обнаружении неисправности происходит перенос данных с тех дисков, на которых замечены сбои. 
Это позволит минимизировать вероятность потери информации. 

При обнаружении угрозы искусственный интеллект анализирует атаку с целью обнаружения её цели. 
Если конечной целью является объект хранилища, то данные с него переносятся на другой накопитель. Таким 


' Легкий способ подделки контрольной суммы с помощью коллизий / хаКер.га. ОВГ: В&рз://хаКер.га /2012/11/22/ 1 4-{аКе-свескзит/ (дата 
обращения: 04.05.2020). 
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образом, в случае успешной атаки, ИИ будет предпринимать эффективные действия для сокрытия информации 
[8]. 

Поскольку с данными работают пользователи, то важную роль играет человеческий фактор. Поэтому 
система устроена таким образом, чтобы с данными могли работать только те пользователи, у которых имеется 
доступ. Здесь есть аналогия с мандатной моделью управления доступом. Исключением является то, что разда- 
чей доступа управляет пользователь, создавший информацию в системе и являющийся её владельцем. 

Заключение. Таким образом, использование искусственного интеллекта при организации хранилища 
данных повышает скорость работы системы. Ограничение доступа пользователей в алгоритме работы системы 
улучшает информационную безопасность. Искусственный интеллект с встроенными алгоритмами машинного 
обучения позволит оперативно реагировать на любую ситуацию, влияющую на состояние системы. Анализ со- 
стояния накопителей позволит избежать возможного аппаратного или программного сбоя. Минимизация чело- 
веческого фактора в работе системы способствует улучшению её функционирования и более глубокому анали- 
зу пользовательских запросов. Кроме того, сбор информации о возможных атаках позволит поддерживать на 
должном уровне безопасность системы. 
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